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INTRODUCTION 


The work accomplished under the subject contract includes 
the following: 

1. Chargeable permanent magnet technology development 

2. Magnetometer development (Mod #4) 

3. Horizon scanner breadboard (Mod #9) 

4. Implementation of a Magnetic Acquisition System (Mod #11) 

5. Development of electronics for dynamic closed loop 
operation of chargeable permanent magnets (Mod #12) 

This report constitutes the final report on this contract and 
includes the Design Study Report on chargeable permanent magnets 
as well as summaries of all other work done under the contract. 
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I DESIGN STUDY REPORT - VARIABLE PERMANENT MAGNETS 
1-. 0 Introduction 

This design study report describes the work done under 
the subject contract on the design of variable permanent 
magnets for use in magnetic balancing and control of earth 
orbiting spacecraft. These magnets can be used instead 
of air core coils or electromagnets in applications where 
the objective is to produce, or eliminate, torque on the 
spacecraft through interaction with the earth's magnetic 
field. Such magnets are particularly useful for on-orbit 
magnetic balancing to minimize the effects of magnetic 
disturbance torques. Once the appropriate moment is 
achieved (by command) , the magnet controller is turned 
off and the moment will be maintained indefinitely with 
no power. 

For magnetic control applications, these magnets can 
have substantial power advantages where larger magnets 
are involved and where frequency of switching is low. 

This report deals with design of the magnet for minimum 
size and weight, and describes the sensing configuration 
used for monitoring the state of the magnet. 
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2 , Q FUNDAMENTALS 


( 2 - 1 ) 

( 2 - 2 ) 


( 2 - 3 ) 


( 2 - 4 ) 

( 2 - 5 ) 


Magnetic Poles 

In treating permanent magnets, it becomes convenient 
to deal with free magnetic poles in analogy to the 
free charges in electrostatics. A unit pole is de- 
fined so that two unit poles with a separation of one 
centimeter in vacuum would experience a force of one 
dyne, and the magnetic field strength (H) is defined 
as the force on a unit pole in this field. The • 
fundamental relations giving the force (F) and the 
magnetic field strength (H) are 

F = m H 

H = m 



where m is the pole strength and r is the distance 
from the pole. When two magnets are brought closely 
together in an end to end orientation, the force be- 
tween the magnets is given by 

F = mi m 2 



provided the separation is small compared to the 
magnet length. 

The effects of magnetic material are measured by the 
intensity of magnetization I — which is the number 
of unit poles per cross sectional area. And because 
the dipole moment is the pole strength times separation, 
the intensity is also the magnetic moment density. 


M = m£ 


I = _M 

volume 

Early experimenters conceived of magnetism as a flow 
of magnetic flux — perhaps because of the patterns 
formed in iron filings by the magnetic field. This 
flow was characterized by lines of induction which 
have the same direction as the magnetic field and 
indicate the flow by their concentration. The induction 
is the sum of the effect due to H and to I and is 
defined by 
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( 2 - 6 ) 


B = H + 4ttI 


where the factor 4 tt occurs because a unit pole gives a 
unit magnetic field at unit distance, and the surface 
of a unit sphere is 4 tt. The effects of material are 
often expressed for linear materials in the more 
familiar expression 

(2-7) B = yH 

2 . 2 Self-Demagnetization and Geometry Effects 

Within a magnetic material, we find an interesting 
effect which is quite important in our study of 
permanent magnets. In figure 2-1, a unit magnetic 
dipole is shown with its lines of induction. As is 
quite obvious, the direction of the magnetic field 
off the end of this dipole is aligned with the dipole 
so that another dipole in this position (shown dotted) 
would find an aiding field. Now consider a position 
to the side of this dipole; the dipole field is 
counter to the sense of the dipole, and a dipole of 
the same sense would find an opposing field in this 
region. The implications are quite clearly that: 
first, magnetic material generates its own demagneti- 
zing field; and second, this effect can be minimized / 
by selecting a long slender geometry. 

The same result may be obtained by considering a 
current loop as the source of magnetic field because 
the field due to a current loop with infinitismal 
dimensions is the same as the field due to a dipole 
of infinitesimal dimensions. The differential mag- 
netizing force due to a differential current element 

(i dll) at a distance r is 

— * 

(2-8 ) dH = i dj. x r 


where r is a unit vector directed from the current 
element to the point of observation. 


The magnetic moment of a current loop may be cal- 
culated from 

(2-9) M = i (Area) 

10 
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It should be realized that the magnetic moment of a 
current loop is defined such that a current loop with 
unit magnetic moment will give unit torque in a field 
with flux density (B) , but a unit dipole is defined to 
give a unit torque with unit magnet field strength (H) . 

Although these two types of moment are dimensionally 
different, they will give the same numerical results 
in air in systems of units where y equals one. 



- \ 

Figure 2>-l Field due to a Magnetic Dipole 
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From this we can readily see that the conversion to 
mks units is one pole cm equals 10 -2 ampere meter 2 . 

Several researchers have investigated self-demagne- 
tization effects, but most magnetic materials are 
highly nonlinear and quantitative information is 
usually obtained from graphed experimental data. 

This data has been summarized to reasonable accuracy 
by assuming that the self-demagnetizing field near 
the center of the magnet is proportional to induction 
and formulating an empirical equation for the propor- 
tionality constant. The proportionality constant is 
taken to be a function of only the geometry and is 
given for long slender bars 1 by 

( 2 - 10 ) 04 = y.49, 

where is the magnet length, D is the width of the 
bar, pD is the perimeter of the cross section, and 
aD2 is the cross sectional area. 

Again, this is an empirical relation giving only the 
self-demagnetizing force near the center of the magnet 
and makes no attempt to include the effect of proper- 
ties of the material on the field distribution, but 
results obtained from this relation have been quite 
accurate . 

Using the assumption that the self-demagnetizing 
force is proportional to induction in the center of 
a magnet, we may determine the induction and effective 
net magnetizing force quite simply by using the load 
line technique. This technique is best described by 
example. 


■^Parker & Studders (Ref, 8) page 165 
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Figure 2-2 Load Line Example 


Turning to Figure 2-2, the procedure is simply to 
draw a load line of slope -(b/H) d through the H axis 
at the applied magnetizing force on a graph of 
B vs H for the material. This line represents the 
locus of all points which satisfy 

B =^B ^ (^applied - H) 


so the intersection of this line with the B-H curve 
of the material indicates both the net magnetizing 
force seen by the material, and the corresponding 
induction. 


B 



Figure 2-3 Typical Sheared B-H Curve 


The same result can also be obtained by redrawing 
the B-H curve of the material for a specific geometry 
so that the self-demagnetizing effects are included 
in the sheared curve. Figure 2-3 shows the effect of 
shearing on a typical curve shape. 


2-5 



( 2 - 11 ) 


( 2 - 12 ) 


Before leaving the subject of self -demagnetization , 
one observation will be noted. This is the effect of 
shearing on the slope of the curve. It is easily 
shown either through differential calculus or plane 
geometry that slopes on the sheared curves can be 
determined from the original slope, and the load line 
slope, from 

cmid (i) D 

which is to say that slopes add like resistances in 
parallel. Very small slopes will be unchanged in the 
sheared curve, but slopes much greater than (B/H)o will 
appear in the sheared curve with a slope of 
(B/H) d . 


(IF) 

' ^ H ^heared 


These considerations will allow us to determine the 
flux density at the center of the bar, but it remains 
to determine the effect of the spatial distribution 
of magnetic intensity (or dipole moment density). As 
our interest here is with macroscopic effects, we need 
only a gross measure as an average. , 

Recalling equation 2-6, if we take averages over the 
entire magnet we have 

Avg(B) = Avg (H) + Avg(4irl) (from 2-6) 

because averaging (arithmetic mean) is a linear 
process. Fortunately, experimental work has been 
done to determine the leakage of induction lines from 
the sides of a bar magnet. From this we can readily 
determine the average flux density, then using the 
load line slope we can determine the self-demagnetizing 
force and if we know the average applied magnetizing 
force we have the average magnetic intensity. 

s 

Figure 2-4 is a graph showing the empirical relation- 
ship between average induction in a bar magnet and 
slenderness ratio X/D, 


2<-6 



AVERAGE INDUCTION 
( % MAX ) 


2 . 


(2-13) 


The significant point is that for J/D ratios greater 
than 4 , the induction averaged over the volume 
of the bar is 70% of the maximum induction at any 
point in the bar. The discrepancy between predicted 
magnetic moment and measured moment for magnets with 
j?/D ratios from 12 to 50 ranged from -1% to +7%. 



Figure 2-4 Effect of Slenderness on 
Average Induction 

Magnetizing Coils 

The standard approximation for the magnetizing force 
due to a coil wound on a long bar is 

H = 4tt ni 

10 17 


where ni is the number of ampere-turns, and J£ c is 
the coil length. 
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c 



In order to determine limits on the validity of this 
approximation, we return to equation 2-8 and calculate 
the magnetic field on the axis of symmetry for several 
current distributions with cylindrical symmetry. The 
coordinate system and applicable dimensions are shown 
in figure 2-5. Note that a uniform radial current 
distribution is assumed for all cases, and that dis- 
tance is measured from the coil center so the coil 
ends are at z = ± J l c 

2 

Y 



Without showing details of the mathematics , the 
results are : 

for a circular current filament 

(2-14) 


H 


= 4£ ni_ 1 +(2z\ 2 
10 D L VD / J 


-3/2 


2^8 





(2-15) 


2 


(2-16) 


(2-17) 


for a current disc (Jl = 0) 



for a cylindrical current sheath (T = 0) 
H = 2jt ni J JL + 2z + X - 2z 

10 I^Jd 2 + (2z+jL) 2 Vd 2 + (2z zO 2 


for a hollow cylinder of current 



The cylindrical current sheath has a particularly 
simple expression for magnetizing force which is 
convenient to establish bounds on the H distribution 
of the hollow cylinder. The field distribution of 
the hollow cylinder will be between the two distri- 
butions given by cylindrical sheathes with diameters 
D and D+2T. 

Inspection of equation 2-16 reveals two useful facts. 
First, the field maximum occurs at the coil center 
(z=0) and approaches 0.4Tmi/£. asymptotically from below 
with increasing length. The length to diameter 
ratio which gives 40 % error is easily shown to be (JJ/ D)=2 
while an (X/D) of 7 gives 1% error. This confirms 
equation 2-13 even for quite short coils. 
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FIELD 


The field at the coil end is very nearly half the 
field at the coil center for as low as 2 and mono- 

tonically approaches half the maximum field as 
(Jt/D) increases. The variation for_ intermediate position 
is shown in figure 2-6. 



Figure 2-6 Magnetic Field Strength Normalized 
to 0.4frni/£c vs position 

While equation 2-17 could be graphed in a similar 
manner, the addition of a third variable (coil 
thickness) provides little additional information. 

4 Magnetic Moment Measurement 

Laboratory measurement of magnetic dipole moment may 
be, in principle, accomplished by either making a 
torque measurement in a calibrated magnetic field, or 
by making magnetic measurements on the magnet. 

The most obvious method is a direct torque measurement, 
but the torque requirements for small satellites have 
historically been in the vicinity of 10 -2 Newton-meters 
(10 -3 in-lb) . It would seem that such delicate torque 
measurements could be circumvented by scaling either 
the magnet dimensions or the ambient magnetic field, 
but both of these approaches have serious drawbacks. 



Because the magnetic core material must be cold worked 
to develop the magnet properties, scaling the core 
dimensions will require the magnetic material manu- 
facturer to scale the size of the cast billet. If 
this is not done, the material must be worked to a 
different degree and the finished stock will exhibit 
different magnetic properties. 

The second approach is even less desirable if a sig- 
nificant scaling is required. Should the ambient 
magnetic field be increased, the magnetic state of 
the core material will be altered. As the core 
material is highly nonlinear, test results could not 
be linearly scaled. 

The magnetic moment of a bar magnet can in theory be 
determined by averaging the flux density within the 
magnet — but this requires detailed knowledge of the 
field distribution within the magnet. 

The flux density within a permanent magnet is not 
readily measured but another method is available for 
determining the average magnetic intensity. 

Consider a magnet of length and the magnetic field 
due to this magnet being measured at a distance R as 
shown in figure 2-7. 


MAGNET 



PROBE 

-X 


Figure 2-7 Position of Magnet & Test Probe 
for Moment Measurement 

If we consider the magnet to be two equal and opposite 
poles of strength m and separation .£ , then the field 
at the measurement point will be 
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( 2 - 18 ) 


( 2 - 19 ) 


m JL 

H= (r 2 +X 2 ) 3/2 

4 

but B = K in air, and the moment is mJl. Therefore, 
we can calculate the dipole moment by making a remote 
measurement of the magnetic field and calculating the 
moment from 



This relation will be in error for any magnet because 
the effective magnetic length is not known precisely, 
and when the measurement is made very close to the 
magnet the effect of local flux density variations 
within the magnet may give significant error. Both 
of these effects diminish rapidly with increased 
measurement distance, with the effect of local varia- 
tions decreasing at least as the fifth power of 
separation . 

Ignoring flux density variations within the magnet 
for the moment, we may determine the sensitivity of 
the measurement to probe placement and to uncertainty 
in the effective length of the magnet. This informa- 
tion can then be used to balance the error caused by 
the magnetic environment against the errors due to 
probe placement and uncertainty in magnet length. 

Consider two attempts at measuring the dipole moment 
of a magnet. Let one measurement be error free and 
the other measurement made with an error 6 in probe 
distance and a relative error y in determining the 
magnet length. The relative measurement error is 
then 


( 2 - 20 ) 


+ ±2 

e =1 4 

(R+6) 2 +( l+Y) 2 j 2 
4 



-1 



2-12 



This may be rewritten to show the position error o 
which will give an error e . 


If 1 \ 2/3 R 2 + JL 2 / 1 V 

'[!+£/ 4 Vl+ £ / 


- (1+Y ) 


Figure 2-8 is a graph of equation 2-21 with zero 
uncertainty in magnet length. This shows the per- 
missible error in probe placement which will cause 
a +5% error in dipole moment when the magnet length 
is known precisely. The significant feature in 
figure 2-8 is the variation of position tolerance 
with separation. If the probe is placed close to the 
magnet to obtain a large field, the tolerance on probe 
position virtually disappears. 



A =30CM 


Figure 2-8 Probe-Position Error to Give 
5% Error in Moment Determination 




Figure 2-9 shows the tolerance on probe position to 
give +5% error in dipole moment when the effective 
magnet length has been estimated 10% too long. This 
amount of error in magnet length will give an error 
of 5% in moment for separation of 1.16 magnet lengths, 
so the position must be exact to obtain 5% accuracy at 
this separation and it is no longer possible to obtain 
5% accuracy for measurement distance less than 1.16 
magnet lengths. 


Figure 2<-9 Probe Position Error to Give 5% 
Error in Moment Determination 
Magnet Length is 10% in error. 
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Quadrupole Effects 


In this section we estimate the effect of non-uniform 
flux distribution on the field strength due to a bar 
magnet. 


Let us approximate a magnet with the array of equal 
charges as shown in Figure 2-10. 




a 


*1 


Figure 2-10 Point Charge Approximation to 
Bar Magnet 

This arrangement of charges may be considered as two 
dipoles superimposed - the dipole formed by the outer 
pair of charges form a dipole of moment itU?, and the 
inner pair of charges form a dipole of opposite sense 
and moment ml?. The net moment is then 2/3 m/. 

3 

The magnetic field due to this arrangement of charges 
is simply the sum of the contributions by the inner 
dipole and the outer dipole. A measurement made as 
in figure 2-7 will indicate a field given by 

mi/3 

(R 2 +JJ2) 3 / 2 
36 

determining the magnetic moment 


m„£ 

(2-22) H = ^ R 2 + £2)3/2 

4 

The relative error in 
from 2-19 is then 
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(2-23) 


e = 1 
2 


R 2 + 1 


12 


R 2 + 1 
P 36 



which is a function of distance relative to the magnet 
length rather than a function of absolute distance. 


Returning to equations 2-22, notice that the variation 
of error with distance will disappear if the inner and 
outer dipoles have the same length. Quite obviously 
the worst variation with distance occurs with a very 
short inner dipole. Taking an infinitismal dipole 
again with half the moment of the outer dipole — the 
error in calculating the dipole moment from equation 
2-19 will be given by 


(2-24) 


£ = ? I 1 "C 1 + ^z ) 37 J 

which is also a function of relative distance. 

Figure 2-11 shows the variation of error with distance 
for these two cases. Both of these assume an inner 
dipole with moment equal to half the net moment. The, 
case of an infinitesimal dipole is the limiting case 
for the relative dipole strengths considered. For 
either case, the error is less than 5% for measurement 
distance greater than two magnet lengths. 


In view of the uncertainty in effective magnet length, 
the selected measurement distance should be as great 
as allowed by equipment sensitivity and by interference 
from the magnetic environment. 

2.6 Characteristics of Magnetic Materials 

The origin of magnetic effects in materials lies in 
atomic dipoles formed by electronic circulation and 
spin. Classical forces tend to align the atomic di- 
poles in opposition to an applied magnetic field, but 
quantum mechanical exchange forces tend to cause the 
dipoles to align in an aiding orientation. The net 
result is that the quantum mechanical forces dominate 
for small groups of atoms which will be aligned aiding, 
but classical forces will limit the size of these 
domains to about 10 -6 cm3 to 10" 2 cm 3 . Macroscopic 
effects then, are due to the actions of these micro- 
scopic domains. 
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The two outstanding features that one notices in 
magnetic materials are explainable in terms of mag- 
netic domains. First is the saturation of the mat- 
erial which occurs when all the atomic dipoles are 
forced into alignment with the field. Because 
magnetic effects are due to the orientation of atomic 
dipoles, the intensity will not increase past the 
point where the magnetizing force causes complete 
alignment of all the atomic dipoles. 

The other outstanding feature is the manner in which 
magnetization changes with applied field. When a 
magnetizing field is first applied to magnetic mat- 
erial, the domains deform elastically. Then as the 
field is increased , the domains in a favorable orien- 
tation begin to grow at the expense of less favorably 
oriented domains, and this process continues until 
either the magnetizing field no longer increases or 
until the material is saturated. This action of 
elastic deformation alternating with domain growth 
gives rise to a multiple valued curve of magnetic 
intensity vs magnetizing field as the field is cycled 
between some minimum and maximum values. 

When a field is applied to a specimen the magnetization 
will begin to increase which will cause the magnetiza- 
tion to increase more rapidly (Figure 2-12A) . As the 
applied field is then reduced, not all of the domains 
will return to their original (random) orientations 
giving rise to residual magnetization. This effect 
is sketched in Figures 2-12A and 2-12B which shows the 
characteristic shape of the hysteresis loop. 

Figure 2-12C shows a minor hysteresis loop which occurs 
when the applied magnetizing field is cycled between 
zero and a value which does not result in domain growth. 
The average slope of this loop is very constant even for 
quite large values of magnetization. This average slope 
on a graph of induction vs applied field (see equation 
6) is termed recoil permeability with symbol yr. The 
recoil permeability can be assumed to decrease about 20% 
with peak magnetization of 90% of saturation. 

Figure 2-12D shows a typical B-H curve shape with 
several recoil loops. Notice the terminology: B m - 
maximum induction; H m - magnetizing force; B r - 
residual induction'; H c - coercive force. 
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From the preceding, we can identify several desirable 
characteristics for permanent magnets: the B-H curve 
of the material should intersect the load line at as 
high an induction as possible in order to minimize 
the mass of the magnet; the maximum required magne- 
tizing force should be low to minimize coil require- 
ments; a long slender magnet should be used to give 
as great a load line slope as possible; and the 
coercive force should be sufficiently high to minimize 
external demagnetizing effects. 
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3.0 THE MAGNETIC CORE 
3.1 Materials 


The primary criteria for selecting the magnetic core 
material are those characteristics that result in a 
minimum weight system. Materials suited to use in 
variable magnetic moment applications are alloys 
with density very near 8 grams/cm3. — and a 10% 
variation is unusual. Thus material density is not 
a significant factor. 

With a choice of materials which have essentially 
the same mass density, the core material which 
gives minimum weight is the one with greatest resi- 
dual induction at the desired operating point. This 
can be shown by combining equations 2-5 and 2-6 to 
give 

(3-25) volume = 4n M 

B-H 

which indicates that the minimum volume magnetic 
core is one with maximum induction. 

The residual induction at the operating point is 
determined by both the B-H characteristics of the / 
material and by self demagnetization due to core 
geometry (see next section) . If we assume that the 
slenderness ratios of concern are greater than 20 
and the residual induction is on the order of 10,000 
gauss, then the self demagnetizing force is in the 
tens of oersteds. This indicates that the optimum 
material has maximum residual induction at a demag- 
netizing force of a few tens of oersteds - but the 
coil requirements should be borne in mind. If a 
material is selected with large residual induction 
but requires a very large magnetizing force, then 
the penalty in coil weight may be very large. 

The core material which gives a minimum weight com- 
bines minimum coercive force and maximum induction 
at a demagnetizing field of a few tens of oersteds. 

Several secondary considerations should be borne in 
mind when evaluating alternative core materials — 
considerations such as available size, special 
treatment required and mechanical properties. 
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Most of the useful core materials must be cold 
worked to develop their magnetic properties, with 
a 95% reduction in area being typical. Because of 
this requirement, some materials are only available 
in thin strip form or small diameter wire. These 
materials also require annealing at approximately 
600°C and must not be deformed after heat treatment. 

A survey of 60 companies was conducted, and the 
three most likely materials for this application 
were selected. Remendur 38, P6 Alloy and Vicalloy 
I were selected for maximum induction over three 
ranges of demagnetizing fields. Remendur 38 has 
the largest residual induction of these three 
materials for demagnetizing fields less than 35 
oersteds, while P6 Alloy has the greatest residual 
induction for demagnetizing fields from 35 to 55 
oersteds, and Vicalloy I extends the range to 100 
oersteds. The properties of these materials are 
summarized in Table 3-1, and the demagnetization 
curves are shown in figure 3-1. 

The choice of material for the magnetic core cannot 
be made without considering the effects of self 
demagnetization, and the requirements placed on the 
magnetizing coil by the choice of core material. 

3.2 Geometry Effects 

As pointed out in 2.2, a permanent magnet has a self- 
demagnetizing field which depends both on the geometry 
and on the level of induction. 

The effect of self-demagnetization is to increase the 
volume of core material required to achieve a given 
monent. This will also increase the coil requirements 
because the magnetizing coil will enclose the volume 
of core material , so a large core will require a large 
magnetizing coil. 

It is convenient to compare core materials graphically 
by superimposing contours of constant residual moment 
per cubed length (M/^3) Q n the B-H curves of the 
materials. This is possible because equation 2-10 
provides a link between the slenderness ratio (C/D) and 
the B-H coordinates, while equations 2-5 and 2-6 provide 
a link between. residual induction, core volume, and 
dipole moment. 

Combining 2-5, 2-6 and 2-10 we obtain 

(3-26) B = 1.024xl0 -4 p2 H 4 

a / M \ 3 

\v) 
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Table 3-1 


Remdender 38 P6 Alloy Vicalloy 




I 


(3-27) 

(3-28) 


On substitution of appropriate geometry factors, this 
becomes 


B = 


1 . 63xl0 -3 H 4 



for a square core cross section 


B 


1.287x10 



for a round core cross section 


Figure 3-2 shows the B-H curves of 3 materials superimposed 
on contours of constant M./Q 3 from equation 3-27. These 
curves facilitate the inclusion of self-demagnetization 
effects when comparing core materials. 


To illustrate the use of these curves, assume that a 
10,000 p-cm magnet is desired with a maximum core length 
of 25.4 cm (10 inches). The ratio M/^ 3 for this case is 
.61, so we follow the M/£ 3 = .6 contour in figure 3-2 
to find that Vicalloy I will require the largest core 
volume with a residual induction of 8100 gauss, P6 alloy 
will require an intermediate volume with 11,900 gauss 
and Remendur 38 will require the smallest volume @ / 
13,600 gauss. The minimum weight core for this 
example would be Remendur 38. 



Figure 3-2 
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3.3 


Interaction with Magnetizing Coil 


The core design will affect the coil requirements 
in two ways: first, through the relationship between 

core and coil geometru; and second, by establishing 
the requirements on magnetizing force. 

The core volume will determine the volume which must 
be supplied with a relatively uniform magnetizing 
force, so the size and shape of the core will obviously 
affect the volume of the magnetizing coil which 
enclose the core. Increasing core length will require 
more wire to cover the coil, while increasing core 
diameter will increase the average length of a turn. 

The magnetizing force which must be supplied by the 
magnetizing coil consists of the magnetizing force 
required to overcome self-demagnetizing effects and 
the net magnetizing force required to saturate the 
magnet. Once the material and required moment are 
specified, the total required magnetizing force will 
be a function of the slenderness ratio. 


The interaction between peak coil power, coil weight 
and magnetizing force is seen by referring to figure 
3-3. This figure shows the shape of the coil weight 
vs. peak coil power curve for a magnetic core of 
fixed size and material. The form of the equation 
plotted in Figure 3-3 is 


coil weight = 


I 2 C 2 

2 

P - P 
H 2 H 2 


C 2 + C 3 


where P is peak power, H is peak magnetizing force 
and the C's are constants (see section 5.2). 


Figure 3-3 shows that there is a threshold coil 
power for a given core design. Below this minimum 
coil power, it is not possible to design a coil which 
will magnetize the core. At large coil power the 
coil weight can be traded one for one with peak 
coil power, but at low coil power the increased coil 
build causes the coil weight to blow up. 

Notice also that" for a given coil design the peak 
power varies with the square of the magnetizing force. 
If the magnetic core is replaced with a material with 
twice the required magnetizing force, then the same 
magnetizing coil may be used but the required peak 
power will increase by a factor of four. 
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Figure 3-3 Coil Weight vs Peak Power 
3.4 Scaling 

Scaling of magnet designs is useful when it is desired 
to produce a different magnet size from a given design 
or when it is desired to evaluate a new design which 
is an inconvenient size. 

It is evident that if a magnet made of homogeneous 
material has its dimensions scaled by a scale factor 
s, then the maximum induction and magnetizing force 
will remain constant. This is consistent with 
equation 2-10 which says that self-demagnetization 
effects depend only on shape and are independent of 
size. 

The spatial variation of magnetic field will have 
the same shape, but measurement distances must be 
scaled . 

This is consistent with the scaling of all distances 
and lineal dimensions. 

The effective pole separation (effective magnet length) 
will be scaled by s, but the effective pole strength 
will be scaled by s . This comes about because the 
pole strength is proportional to the cross sectional 
area of the magnet rather than to the volume of 
magnetic material. 

Finally, the volume and weight of the material, the 
dipole moment, and the energy required to magnetize 
the material will all be scaled by s^ — which is 
what would be intuitively expected. 


3-6 



Summarizing these scale factors, we find that: 
all dimensions and distances are scaled by s;volume 
weight, and moment are scaled by s ; and pole 
strength is scaled by s . 
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GENERAL ELECTRIC PERMANENT MAGNETS • THERMISTORS • THYRITE® VARISTORS 


DESCRIPTION: 

P-6 Alloy combines the advantages of a high hysteresis loss and low magnetizing force into 
one magnetic material. These two important features plus high residual induction and high 
permeability make P-6 Alloy especially useful for hysteresis motor applications. 

This material is supplied in strip or wire form. The strip material is available in widths up to 
4.0 inches and thicknesses from .005" to .100". Wire may be obtained in diameters ranging 
from .005" to .100". Magnetic certification is furnished with each shipment of material. 

CHARACTERISTICS: 

A. Composition: 


Nickel 

6.0% 

Vanadium 

4.5% 

Iron 

44.5% 

Cobalt 

45.0% 


B. Magnetic Properties *(D.C.): 


Maximum permeability 285 

Peak magnetizing force, oersteds 250 

Maximum induction, gausses 16,000 

Residual induction, gausses 14,000 

Maximum energy product (B^H^), gauss oersteds 550,000 

Coercive force, oersteds 60 

Hysteresis loss, joules/invcycle .37 

peak H = 180 ampere turns/inch ^-s.^5'5' 

* At 25°C after age-hardening 
C. Physical Properties: 

Hardness prior to age-hardening 40 R.C. 

Hardness after age-hardening 65 R.C. 

Specific gravity 8.16 

Density, Ibs./cu. in. 0.285 

Electrical resistivity, microhms-cm 30 

Coefficient of expansion, in./in./°C 11.1 x 10 6 

HEAT TREATMENT: 


After stamping or coiling to the desired form, the material should be either vapor degreased or 
emulsion cleaned prior to final age-hardening heat treatment. The material should then be 
placed in an electric furnace with hydrogen or cracked gas atmospheres and held at 600° C for two 
hours. This final heat treatment is necessary to obtain the desired magnetic characteristics of 
the material. 

DESIGN CONSIDERATIONS: 

A. General: 

P-6 Alloy is capable of supplying superior hysteresis characteristics in designs where 
a maximum of 100 oersteds magnetizing force is available. 

The high aging temperature of P-6 alloy allows it to be subjected to a much higher 
operating temperature than the chrome, tungsten and cobalt magnet steels without 
resulting in subsequent aging. The material is magnetostrictive, however, and lamina- 
tions should therefore be stacked and contained with a minimum of stress applied to the 
material. Figures 1, 2 and 3 respectively show the hysteresis loop, energy product and 
magnetization curves; the hysteresis loss vs. maximum flux density; and the hysteresis 
loss vs. magnetizing force for fully aged P-6 Alloy. 
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4.0 


MAGNETIZING COIL WEIGHT 


In this section we treat the problem of determining the 
wire weight for a magnetizing coil. 

4 . 1 Derivation of Coil Weight Equations 


The resistance of a length of wire in terms of its 
resistivity and dimensions is given by 


(4-29) 


(4-30) 


(4-31) 



7rd^ 

or the resistance for a specific size and material 
is given by 

R = ri w 

where r is lineal resistivity - e.g. resistance per 
thousand feet of wire length. Similarly, the weight 
of a length of wire in terms of wire dimensions and 
density is given by 

W w = JL w Udi <$w 

4 

or an expression similar to 4-30 can be written for 
wire weight. 



Figure 4-1 Coil Dimensions 



(4-32) 


(4-33' 


(4-34) 


The length of wire required to wind a coil as shown 
in Figure 4-1 can be approximated in three steps. 

The first order approximation is that wire length 
equals the product of the perimeter of the cross 
section and the number of turns. 

-£ w = (pD) n 

The second order approximation corrects for the 
effect of coil build - i.e., for multilayer con- 
struction. Because alternate coil layers will have 
opposite pitch, the diameter of the coil will increase 
by two wire diameters for each additional layer of 
wire. The length of one turn in the jth layer will 
then be 

(pD) + 2Trd (j - 1) 

2 

and the length of all turns in X layers is approxi- 
mately 

j£ w = n (P°) fl + 1 xd *7 

L p D - 1 

It is of value to note that (u/p) is nearly 1, so the 
second term in the brackets is nearly the ratio 
of coil build to coil diameter. 

The next refinement is very rarely needed unless the 
wire diameter approaches the core diameter. The 
additional wire length due to the pitch of the winding 
is simply the length of the coil multiplied by the 
number of coil layers. The expression for the length 
of wire in a coil of X layers is 

& w = n (p°) r 1 + i ^.1 + ^^c 

p D J 

This has been found to give ±1% accuracy for care- 
fully wound coils. From this we can 

compare the amount of wire used for different cross- 
sectional shapes. For instance, consider two cores 
with equal cross sectional area. One having square 
cross section with side w, and the other being 
circular with diameter D. 
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For the square : 
£wl = 4nw + mrXd 


For the circle : 


Xw2 = nD + n-irXd 


Invoking the equal area constraint we have for the 
ratio of the wire lengths 


Iwl 

iw'2 



This ratio is plotted vs (A/d^) with \ as a parameter 
in the graph of Figure 4-2. This shows that the 
penalty for using a core with square core cross 
section rather than a round cross section varies from 
13% when a single layer of fine wire is wound on a 
large core to 26% when many layers of heavy wire are 
wound on a small core. / 

Included in the figure is a nomograph to convert cross 
sectional area and wire gauge to (A/d 2 ) , The nomo- 
graph assumes double formvar insulation, so the 
diameter used when constructing the nomograph was 
greater than the bare wire diameter listed in wire 
tables. 

The use of the nomograph is illustrated by taking a 
core with 0.5 square inch cross sectional area with 
10 layers of 24 AWG. The intersection of the 24 AWG 
line and the 0.5 in 2 line gives an (A/d 2 ) of 1005. 
Following the (A/d 2 ) line down to the 10 layer curve, 
we find the length ratio is 1.160. For this example 
16% more wire is needed for a square cross section 
than for a circular cross section of equal area. 

Figure 4-3 facilitates comparison of wire length with 
circular cross sections when only the wire diameter 
and the number of layers is changed. If the same 
core is used, then the periphery (pD) will not change 
and if the number of turns is also constant, then the 
relative wire length is 



RELATIVE WIRE LENGTH 


I 


Fig . 4-2 Wire length for Square Core Relative to Wire length 
for a Circular Core 

J-w_ = 1 + Ad 

npD D 

from equation 4-33. 


To illustrate the use of Figure 4-3, assume that a 
coil is to be designed with 1000 turns of copper wire 
and 5ft resistance on an 8 inch long by .5 inch diameter 
round core. 






WIRE TABLE 


AWG 

Diameter 

in Inches 

Ohms per 1000 

feet 


Bare 

Double 

Formvar 

Aluminum 

Copper 

10 

.1019 

.1055 

1.637 

.9989 

12 

.0808 

.0842 

2.603 

1.588 

14 

.0641 

.0673 

4.139 

2.525 

16 

.0508 

.0538 

6.584 

4.016 

18 

.0403 

.0431 

10.47 

6.385 

20 

.0320 

.0346 

16.64 

10.15 

22 

.0253 

.0277 

26.46 

16.14 

24 

.0201 

.0223 

42.08 

25.67 • 

26 

.0159 

.0179 

66.90 

40.81 

28 

.0126 

.0145 

106.4 

64.90 

30 

.0100 

.0116 

169.2 

103.2 

32 

.0080 

.0094 

269.0 

164.1 


TABLE 4-1 


The periphery of the core is 1.57 inches (.5tt), 
so 1000 turns will require at least 1570 inches 
or 131 feet. 

The lineal resistivity is then on the order of 5ft 
per 131 feet or 38ft per 1000 feet. As we expect 
to use more than 131 feet of wire because of coil 
build, we look up the wire size with a slightly 
lower lineal resistivity in Table 4-1 and find 24 
AWG with 25.6712 per 1000 feet. 

Next we approximate the number of layers "by taking 
the length of a single layer coil (number of turns 
multiplied by wire diameter) and divide by the 
desired coil length. 

\ i 
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1000 turns x .0223 inch/turn 
A - 8 inch 

A * 2.8 

With the preliminaries completed, we use the wire 
nomograph at the top of figure 4-3 to find that (D/d) 
for a .5 inch diameter core and #24 AWG is 22.4. 

Next, enter figure 4-3 at (D/d) = 22.4 and 3 layers. 
The wire length ratio is 1.13, so the wire length is 

j0 w = 1.13 x 131 feet 

iw - 147 feet 

and the resistance of the coil is 
Rc = 147 feet x 25.67^/1000 feet 
Rc = 3.77ft 

If this coil is wound with full layers rather than 
allowing the last layer to partially cover the coil, 
then the coil length will be 

Jic = 1000 turns x .0223 inch/turn 

3 layers 

X c = 7.43 inch/layer 

Following this procedure we can quickly evaluate 
alternative wire sizes, wire materials and determine 
the fraction of the coil length occupied by the last 
layer . 

From Figure 4—2 , a considerable savings in wire length 
results by replacing a square core with one of cir- 
cular cross section. It is inconvenient to fabricate 
a magnet with a true circular cross section from 
stacked laminations, but a circular cross section may 
be approximated by varying the widths of the lamina- 
tions . 

If many laminations are used, the cross sectional 
shape may be approximated by a regular polygon of 
(2N+2) sides - where N is the number of laminations. 
This accounts for 2N exposed edges of the laminations 
and the two exposed sides. 
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This polygonal approximation implies a variation 
in the thickness of the laminations because the 
sides of the polygon are not perpendicular to the 
width of the laminations. The ratio of the minimum 
thickness to the maximum thickness is 


«£sin 71 
N+l 

T = minimum lamination thickness 


T = maximum lamination thickness 
Ji = length of a side of the polygon 
N = number of laminations 


With this limitation in mind we proceed as follows. 

The perimeter and area of the polygon are given by 

Area = (2N+2) D 2 tan tt 

4 2N+2 


Perimeter = (2N+2) D tan ir 

2N+2 



D = diameter of inscribed circle 

N = number of laminations 

2N+2 = number of sides to polygon 

But if N is greater than 5 , then 

(2N+2) tan/ ir \= 12 tan (15°) 

12N+2/ 


= 1.023tt 

So, within 2%, the area and perimeter of the 
polygon are the same as that of the inscribed circle. 

From this we conclude that the advantage in wire 
length for circular cross section can be gained by 
stacking as few as five laminations if the widths are 
cut to approximate a circular cross section. 
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Conductor Materials 


At present, the only practical conductor materials are 
aluminum and copper , but should another material be 
considered, some of the most desirable properties are 
given here. 

The conductor characteristics of prime importance are 
density and resistivity. Quite obviously, the conductor 
material should be selected for minimum density to 
minimize coil weight, and should be selected for minimum 
resistivity to minimize power. Unfortunately, these 
two characteristics do not necessarily go hand in hand, 
and the optimum material will depend on the specific 
application. In the case of a planar' air coil with 
little build, the power and weight may be traded one 
for one. In this case the product of density and 
resistivity becomes a meaningful measure of conductor 
quality. But when the coil build is appreciable, this 
is no longer a useful guide. 

It is possible to predict the crossover point where 
comparison of the resistivity-density product of two 
materials will no longer indicate the conductor material 
for a minimum weight coil , but this must be done for a 
specific core geometry and specific core material. 

(See section 4.3) 

Should a conductor material other than aluminum or 
copper be considered, several other properties besides 
resistivity and density must be considered. Copper 
and aluminum have well known properties and any other 
material will require extensive effort to become 
established. The properties which must be established 
include flexibility, mechanical stability under winding 
pressure and tension, susceptibility to creep, effects 
of impurities, and other odd properties such as the 
tendency of silver to migrate through dielectric 
materials or the growth of crystals from the base 
metal (whiskers). 



Coil Power 


The magnetizing coil for a given magnetic core may be 
designed for a wide range of peak magnetizing power, 
with coil weight decreasing as peak power is increased. 
The lower limit on coil power is fixed by the magnetic 
core design and the coil conductor material. These 
effects are easily demonstrated by taking specific 
examples. 

The minimum power effect is illustrated by examining 
the effect of adding more layers of wire to the magnet 
assembly in Figure 4-4 which is a single layer coil 
wound with wire whose diameter is one tenth the core 
diameter. The change in coil resistance due to the 
additional coil layers is proportional to the increase 
in wire length . 

The original layer of wire will have a mean length per 
turn of (d + lOd) , and the mean turn length in a second 
layer will7Td(10 + 3). The mean diameter of subsequent 
layers will increase by two wire diameters for each 
additional layer so the mean length of a turn in the 
n th layer is1Td(9 + 2 A ). / 
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Using 

P = i 2 R 

and letting subscripts denote the total number of 
layers, the power in the original coil is given by: 


Pl = (ij .) 2 Ri 

Now when we add a second layer of wire we will have 
doubled the number of turns so we must halve the 
current to hold the magnetizing field strength 
constant. 


P 2 = ii 2 R2 
2 ~ 


= 2 R x 1 + Trd (10+3) 

2 Trd (10+1) 


P 2 = (i x ) 2 Ri x 1.09 
2 

For a three layer coil we have 

P 3 = ii 2 r 3 
3 

= ±i 2 R x 1 + ud (10+3) + ird (10+5) 
3“ rrd (10+1) ird (10+1) 


P 3 = (ii ) 2 R^ x 1 . 18 
~3 

and for a X layer coil the power is 

PX = (ii) 2 Ri 10+X 
X 10+1 
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Now we have obviously reduced the power requirements, 
but not as rapidly as we have increased the total 
number of turns. The limiting power with an infinite 
number of coil turns (and layers) will be given by 

Lim PA = (ii)^Ri 1 

U /d +1 

A-*- 00 


This would seem to indicate that there may be an 
advantage to using fine wire so that the limiting 
power would be less. But the actual effect is only 
to increase the power for a single layer coil so 
that a greater reduction in power is possible. 

Notice in the previous illustration that coil 
weight increases more rapidly than the number of 
turns because of the increased length of a turn. 

By the same line of reasoning used in determining 
the coil resistance, the coil weight is 

WA = W-i A 10+x 

10+1 

I 

and the product of coil weight and power becomes 

PAW A = PjWi 10+A 2 
10+1 

or to give specific examples 
P 2 W 2 = 1.19 P]Wi 


P 3 W 3 =1.40 PjWi 


P 4 W 4 = 1.62 Pj^ 
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Figure 4-5 shows the effect of increasing the number 
of coil layers on the coil weight, coil power and on 
the coil weight - power product. These results are not 
restricted to this specific geometry, but can be gener- 
alized to show that the minimum magnetizing power is 
independent of the wire diameter. 



Layers (X) 


Figure 4-5 Variation of Coil Parameters with Coil Layers 
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In the limit when extremely large coils are used, the 
actual shape of the core cross section becomes im- 
material, and the limiting power is a function of the 
magnetic field strength required, the length of the coil 
and the resistivity of the conductor. By combining 
equations 2-13, 4-29, 4-31 and 4-33 we may solve for the 
limiting power with an infinite amount of conductor. The 
limiting power is given by 


P Q = 2.81 H pk 2 £c r 

which for aluminum and copper conductor materials 
becomes 


ai 

7.95 

x 10“ 6 

H 2 
pk 

slc 

cu 

4.85 

x 10" 6 

H . 2 
pk 

SL C 


where a reasonable estimate of insulation thicknenn 
has been included. 

These results point to another effect due to coil 
build. When a coil is constructed of aluminum and 
another coil is constructed of copper, with the same 
wire diameter, number of layers and turns, the alumi- 
num coil has one third the weight because aluminum is 
less dense. The copper coil requires 60% as much power 
as the aluminum coil because copper has a greater con- 
ductivity. We might therefore expect to be able to 
construct an aluminum coil that requires the same power 
as the copper but with less weight. But coil weight 
and power do not trade one for one, so an aluminum coil 
is lighter only when the power is above some threshold 
value. 


This threshold power can be determined by setting the 
weight of two coils equal and solving for the power. 
The resulting threshold power for materials denoted by 
subscripts 1 and 2 is approximately 


P, . = 2.81 x H . ' 
1-2 pk 


£c 


_L - JL 

IT 12 

1 _ 1 ~ 

p 16 1 p262 


4-14 



And this threshold power for the most common choice 
of aluminum versus copper conductor is 


P . =26 

al-cu 


x 10 -6 H , 2 £c 
pk 


For a 10,000 p-cm magnet 25 cm long and requiring 150 
Oersteds, this threshold power is on the order of 
fourteen watts. 


The peak coil power for a permanent magnet is useful 
for determining the weight of the magnetizing coil and 
for determining the circuit requirements for the coil 
driver. However the peak coil power is not a valid 
basis for comparing the power of different systems un- 
less the time duration of the power drain is considered. 


The charging energy — or average power during a charg- 
ing cycle can be shown to be independent of the coil 
resistance. This can be seen by considering a series 
RL circuit (Fig. 4-6) . When the switch is closed the 
current willfollow the familiar exponential rise as 
shown in Fig. 4-6. 

If the circuit is interrupted when the current reaches 
some portion of the maximum possible V q , 

R~ 




then the energy supplied will be the product of the 
supply voltage and the average current supplied to the 
circuit. If the resistance is then doubled, the supply 
voltage must be doubled to maintain the same maximum 
current. Now, when the switch is closed the current 
will rise twice as fast because the L/R time constant 
has been halved. If the switch is opened at the same 
current level as before, then the peak current will be 
the same as before, but the energy supplied to the 
circuit has not changed. 
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This same line of reasoning can be extended to show that 
the energy required to apply a given magnetic field to a 
given core is independent of the coil parameters. This, 
of course, assumes that the spatial distribution of the 
magnetic field is the same for all coils. A very large 
diameter coil will obviously have more leakage than a 
small diameter coil of the same length, but when treating 
relatively small coils the energy required to charge a 
magnet is independent of coil turns, wire size and con- 
ductor resistivity. 


4.5 SCALING 


When the magnetic core is scaled (see 3.5) to produce 
a different magnet size from a given design, the coil 
dimensions must also be scaled. In this section we take 
a scale factor s, scale all coil dimensions, and determine 
the effects of this scaling. 

If the coil length and wire diameter are scaled by the 
same factor, then the number of turns which will fit on 
one layer will be constant. Similarly, the number of 
layers is constant when the wire diameter and coil built 
are scaled by the same scale factor. 


The length of a turn in any particular layer is scaled by 
s, and consequently the wire length is scaled by s. 


The resistance of the coil varies directly with the wire 
length and inversely with the square of the wire diameter, 
so the coil resistance is scaled by s~l. 


For the magnetic field strength to remain constant, the 
current must be scaled by s. This can be seen from 
equation 2-13. 


H = ... 4irni.ii = . 4.ira2-i2 


*1 


*2 


..Ami-ii - . 4irnii2 
£l S 2-1 

i 2 = si! 
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The scale factor for peak power can now be readily deter- 
mined from 

P = i 2 R 

by using the scale factors for current and resistance. 

The scale factor for peak power is s. The energy required 
to charge the magnet is proportional to the volume of the 
coil and core - e.g. to s3. Thus peak power scales by s 
and magnetizing time by s 2 . 

The last scale factor to consider is for voltage. This 
is easily found by taking the product of the scale factors 
for current and resistance/ with the result that voltage 
is unchanged. 

Another simple scaling method is possible - the coil 
dimensions may be scaled by s , but the wire diameter is un- 
sealed. This method of scaling gives a scale factor of s 2 
for coil resistance and wire length. The scale factor for 
peak power, number of turns per layer and number of layers 
is s. The coil voltage is scaled by s 2 and peak current 
by s“l. The magnetizing time is scaled by s 2 , and peak 
power by s so that magnetizing energy is again scaled by 
s3. 

/ 

* ] 

There are two limitations to coil scaling: fine wire 

usually has relatively heavy insulation; wire diameter is 
quantized at approximately 20% intervals. Neither of 
these is a serious objection and only a minor adjustment to 
a scaled coil design is ever needed in practice. 



5.0 MAGNET WEIGHT 


5. 
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Core Weight 

In section 2.2, it was pointed out that short magnets 
generate a self-demagnetizing field which partially 
demagnetizes the magnet in the absence of an externally 
applied magnetizing force. This geometry effect was 
used in section 3.2 to find the maximum residual 
induction in the magnetic core as a function of the 
moment and core length. 

Because the dipole moment, induction and volume of 
core material are related by equations 2-5 and 2-6, 
and by assuming 30% leakage, we may write for the 
core volume 

Volume = 4 M r 

UTTB' r 

where M r and B r are the residual moment and the 
residual induction respectively. The weight of the 
magnetic core is simply the product of core volume 
and the density of the core material, so data points 
from 3-2 may be used to generate curves of core [ 

weight vs residual moment vs core length. 

Figure 5-1 gives core weight as a function of core 
length for a 1000 pole-cm magnet. As might have been 
expected, these curves flatten out for long magnets 
where the self-demagnetizing force becomes small. 

This curve may be used for other magnet sizes by 
scaling (see section 3.5 or 5.4). The weight will 
scale directly with the moment and the length with 
the cube root of the moment. To illustrate, consider 
the problem of finding the minimum core weight for a 
12,000 pole-cm P6 magnet in 25cm of length. The 
dipole moment is scaled by a factor of twelve, so the 
length is scaled by the cube root of twelve or 2.29. 

The length scaled down to 1000 pole-cm is 10.9cm 
(25cm divided by 2.29,) so we enter figure 5-1 at 
10.9cm. The indicated weight is 12.9 grams for P6 
material, which must be scaled by twelve. The minimum 
weight for a 12,000 pole-cm magnet of P6 material in 
25cm of length is 155 grams. 
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Coil Weight 

The weight of the magnetizing coil is determined by the 
magnetic core material, the required dipole moment, 
magnet length and the peak power . 

Specification of the magnetic core material, dipole 
moment and magnet length will determine the slenderness 
ratio and the core volume. These in turn determine the 
peak magnetizing force necessary to satisfy the material 
requirements and to overcome the self-demagnetizing 
force (section 2.2). The coil weight is an increasing 
function of core volume and peak magnetizing force, 
but is a decreasing function of coil power. 

The peak coil power cannot be reduced beyond a 
threshold power because coil losses increase rapidly 
with coil size (section 4.3). The minimum coil power 
is given by 

P Q = 2.81 H pk 2 /c 


Several expressions for the coil weight may be derived 
from the information in chapters 3 and 4 . One convenient 
form giving the coil weight relative to the core weight 


is 

W wire 

W 

core 


ff 2 c 
a 

4.437 P 



^wire 

^core 

8£ + 4.452 

9 



This equation with the numerator set equal to one 
is plotted in Figure 5-2. Notice that for large 



1 2 5 10 20 100 



Figure 


5-2 Variation of Coil Weight with Peak Power 


peak power the curve is asymptotic to (4.437 'P/P ) 


-1 
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5.3 


Magnet Assembly Weight 


The weight of the total magnet assembly consists not 
only of the weight of the core and coil, but also 
includes the weight of the housing, potting or 
encapsulation, connector (s) , and any internally, 
mounted sensors for temperature, moment, etc. 

The sensors will generally be very small, light 
devices and their weight will be insignificant. But the 
weight of the housing, potting and connector will 
depend on the specific circumstances and can be 
expected to increase the total weight by at least 25%. 

In view of the number of variables, it is more 
convenient to find the weight of the core and the 
weight of the foil separately as detailed in 
sections 5.1 and 5.2. 

5 . 4 Scaling 


When scaling an extant magnet design to obtain a different 
size, there are two principle limitations. The most 
troublesome of these is that wire is available only 
in discrete sizes, so exact coil scaling will not ! 
usually be possible. 

The second important limitation is the effect of scaling 
on the magnetic properties of the core materials. As 
discussed in 3.5, this is not a serious limitation if 
the magnet core is fabricated of stacked limitations. 

Table 5-1 lists most parameters of interest with their 
scale factors for a magnet with all dimensions scaled 
by a factor s. These scale factors are exact for 
materials whose properties do not change with size - 
in spite of the heuristic development in sections 3.5 
and 4.5. 
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TABLE 5-1 MAGNET ASSEMBLY SCALE FACTORS 


Fixed 

Voltage 

Scaling 


Fixed 
Wire Size 
Scaling 


Moment (M) 

s3 

Total Weight 

s3 

Magnetizing Energy 

s3 


s-3 

s3 

<=3 


Core length ( ) s 

Core diameter (D) s 

Distance for equal remote s 

field (r) 

Magnetic field strength (H) 1 

Magnetic induction (B) 1 

core weight s^ 


s 

s 

s 


1 



coil length ( c) 
coil turns ( ) 
coil layers ( ) 
turns per layer 
wire diameter 
wire length 
coil weight 


s 

1 

1 

1 

s 

s 

s3 


s' 

s 

s 

1. 


coil current 
coil resistance 
voltage 
power 



s 

s 

s 

s 


-1 

3 

2 


where s is a dimensionless scale factor. 
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6.0 Magnet Design Curves 

The following curves can be utilized to estimate the size, 
weight, and power of chargeable permanent magnets. The 
design curves are based on the elliptical geometry described 
in the next section. 

6. 1 Core Weight and Length 

The size and weight of the magnetic core for a wide range 
of dipole moments can be determined from Figures 1 and 2. 

6 . 2 Coil Weight 

The weight of the magnet coil will depend on the peak 
power used to change the magnetic state. By peak power 
is meant the power required to just saturate the magnet. 

This dependency can be seen by performing a mental 
experiment. If half the turns are stripped off the coil, 
the resistance, wire length, and turns will be approximately 
halved. If this stripped coil is connected to the same 
DC voltage as used by the full coil, the current will 
double but the nimber of ampere-turns will remain constant. 
This stripped coil will thus have approximately the same 
effect on the magnetic core as the full coil but at twice 
the peak power and half the weight. 

Figure 3 shows the tradeoff between coil weight and peak 
power for several magnet sizes. These curves do not 
include the effect of insulation on coil weight. 

Notice that for large peak power, the coil weight and 
powercan be traded one for one (i.e., doubling the coil 
weight reqires half the power) , but small peak power 
results in a considerable weight penalty because of coil 
build-up effects. 

Insulation weight will range from 5% of coil weight for 
high power coils to 15% for low power coils. 
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Figure 1. Length of Core vs. Maximum Dipole Moment 
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Figure 2. Core Weight vs . Moment 
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6 . 3 Switching Energy 


It is somewhat misleading to talk of power used by 
variable permanent magnets, as power is consumed only 
when the dipole moment is changed. A more meaningful 
concept is the energy used for each change. This 
energy goes to three sinks in the circuit, but must all 
come from the satellite power supply. 

The first energy sink is the hysteresis of the magnet 
material. Any change in magnetic state requires energy 
proportional to the area of the B-H loop traversed. For 
well designed systems, this is the dominant energy sink. 

The second sink is I 2 R loss in the coil. 

The final sink is the charging circuit. The simplest and 
most efficient charging technique is to switch the coil 
directly to the spacecraft power supply line, in which 
case this loss is essentially zero. For most magnets, 
assuming there are no severe limits on the peak current 
available from the power supply, this is by far the 
preferred technique. It is possible to isolate the space- 
craft power supply from these current transients by using 
either capacitors or rechargeable batteries , but this 
carries a definite system weight penalty. 

Figure 4 can be used to estimate both the time and the total 
energy required to switch a magnet between its extreme 
states. This includes core and coil losses. From 
these curves it is evident that approximately 6 x 10 5 
watt sec/pole cm is required to change the state of a 
magnet, more or less independent of magnet size. Thus 
a 10,000p-cm magnet can be switched between its extremes 
(M = 20,000 p-cm) at a cost of 1.2 watt sec of energy. 

6.4 Example 

The following example illustrates how to use these curves 
to determine the approximate magnet size and weight for a 
given set of requirements. 

Maximum dipole moment required: 12,000 p-cm (12 amp m 2 ) 

Power supply voltage: 25 V 

Peak currnet allowable: .3 amps 

From Figure 2, the weight of the core will be .2 lb. 

The length of the magnet assembly as taken from Figure 3 
is 15 inches. 
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To determine the coil weight, we use Figure 4 and estimate 
the position of the 12,000 p-cm curve. We estimate the 
coil weight for 7.5 watts and 12,000 p-cm as 0.4 lb. 

The magnetizing time for the magnet can be estimated from 
Figure 5. For a 12,000 p-cm magnet at 7.5 watts this 
time is approximately .14 seconds. 


The estimated parameters are tabulated below for comparison 
with an actual magnet assembly. 

Estimated Actual 


Maximum Dipole Moment 
Core weight 
Coil length 
Coil weight 


12,000 p-cm 
.2 lb 
15 inch 
0.4 lb 


12,300 p-cm 
.22 lb. 

15 inch 
0.48 lb 



COIL WEIGHT(lbs) 
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7.0 Description of Final Magnet Design and Performance 

7.1 Introduction 

This section presents test results and construction 
details for an elliptical chargeable permanent magnet. 

The test magnet is in the form of a prolate ellipsoid 
with integral magnetizing coil. 

7.2 Sensing 

One of the most difficult design problems is the imple- 
mentation of the sensor for determining the state of the 
magnet. Such a sensor is desirable in balancing 
applications for telemetry purposes, and is essential 
in control applications because there is no way to 
properly control the state of the magnet without feedback 
from such a sensor. 

The most convenient sensing device is the Hall sensor, 
which may be attached directly to the magnet. The primary 
difficulty arises from the fact that the magnetic state in 
the vicinity of the sensor does not necessarily represent 
that of the entire magnet. The accuracy of the sensed 
moment is therefore dependent on magnet history. The 
optimum configuration of the many that were investigated 
was for the sensor to be mounted near the end of an 
elliptical core, with a soft iron pole piece outboard of 
the sensor. This configuration is described in the 
drawings, and gives the sensing accuracy described in the 
accompanying curves . 

7.3 Test Results 

The magnet core is a 14 inch long by 0.34 inch diameter 
prolate ellipsoid and the weight of the magnet core and 
magnetizing coil is 0.575 pounds (262 grams). 

The maximum residual magnetism for this magnet is 13,400 
pole-cm and accuracy of the telemetry is *200 pole-cm 
(±1. 5%) . The accuracy increases when the magnet is 
utilized to less than its maximum moment capability - 
but a large decrease in moment is accompanied by a small 
increase in accuracy. These effects can be seen in the 
sensed moment vs actual moment curves of figures 1 and 2. 


7-1 



5s::ss:s:::sK:s:::s:s:::ssssh:rir.rTM:2Si , iSrr^^^PM;ri^er^i^^£35SH:ss::ss: 

SSSaBSUSSSSSaBSgSaSBSS^BmRraSBBmRB^^SBraSSSSSSSBBS 

I i 'g[ i i ' Hinillllli ■■■■■HBHBBU 

■sisK^IES^si^S 

HHfi H HHHHHHH HHBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB*SBBBBBBBBBB 
SSSSSSS8SSSBSBBSSSE3SSBSSSSBSSSBSSSBS8SSSSSSSSSSSSB3SSSSSSSSSS SBSSSSSSSS 
SSSSSSBBSSSSSSBSSBBBSBBSBSSBBSSBBSSSSSSSBSSBSBSBS:BBSSBSSBSS8SS.BSBBSSB&SBS 


■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■I 

SSSSSBSSSBSSSSSSSSBBSSSBSBSBSSS: EBSSS 


■■■■■■■ ■■■■■ I 

ISSnSSnSil 

■■■■■■nnl 



I ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■»■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■! 
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■I 


2 




7 . 4 Magnet Construction 

The test magnet was fabricated in the same major steps 
as are all ITHACO magnets. These steps are: fabricating 

the magnetic core which includes bonding laminations into 
a stack and grinding to final form; magnet cap assembly; 
winding the coil ; potting and housing the complete 
magnet. 

7.4.1 Magnet Core 

Construction of magnetic cores begins by cutting 
Remendur 38 core material to rough size from 0.02 inch 
strip stock. The rough cut strips are heat treated at 
1100°F in an inert atmosphere. After a two hour bake 
the strips are cooled at 100 °F per hour to 900 °F and 
then removed from the oven to cool in ambient air. 

The heat treated strips are next bonded together 
with epoxy to form a rectangular bar and then ground to 
final shape. This magnet was ground to an elliptical 
shape with one end cut off square as shown in figure 3a. 

7.4.2 Magnet Cap 

The magnet caps are 3 piece assemblies as shown 
in sketch SK00224. Part #1 in this sketch is used as a 
partial coil form on both ends of the magnet. Figure 3b 
shows parts #1 bonded to the magnet core. 

One end of the magnet core is ground flat to provide 
a flat face for the Hall sensor. The sensor is a small 
rectangular plate and is bonded to SK00224 part #3. Part 
#3 is then slid into part #2 with the sensor facing out, 
the sensor leads are placed in the wire channel and this 
assembly is slid into part #1. The sensor is then 
sandwiched between part #3 and the flat end of the magnet 
core. The final assembly step for the magnet cap is 
securing the plug part #2 in place with set screws. The 
end caps and core are now ready for winding the coil. 

7.4.3 Magnetizing Coil 

This magnetizing coil is a conventional solenoid 
consisting of 3703 turns of #26 AWG copper wire in 5 layers. 
After winding, the coil was coated with PC22 polyurethane 
and cured. The magnet assembly is shown in figure 3c before 
conformal coating. 
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At this point the magnet assembly is functionally 
complete. The remaining operations are inserting the 
magnet into its housing, potting the assembly and painting 
the housing. 

7.4.4 Potting and Housing 

The housing selected for the test magnet was a 
15 inch length of fiberglass tubing with an inside 
diameter of 5/8 inch. The space between the magnetizing 
coil and the tube was filled with PR1524 polyurethane. 

This complete magnet assembly is 15 inches long 
x 3/4 inch diameter and weighs 1.0 pounds including a 
6 ft. cable. 

7.5 Test Method 

The magnet was tested by placing the magnet on a wood 
table and sensing the remote field of the magnet with a 
flux gate magnetometer. The magnetometer probe was 
positioned parallel to the long axis of the magnet at 
a spacing of 45.5cm. At this spacing the calculated 
scale factor is 10 milligauss per 1000 pole-cm. 

The magnetizing coil was energized for this test by a 1400uf 
bank of capacitors in the MMCA Test Panel. The MMCA Test 
Panel also has a polarity reversing switch and a selector 
switch which connects the capacitor bank to either an 
external power supply or to the magnetizing coil. 

An amplifier was used to increase the output voltage of 
the Hall sensor to a more convenient level. As shown 
in figure 4, this consists of a differential amplifier 
with a gain of 50, which provides an offset adjustment 
control and also includes a bias network for the Hall 
sensor. 

The entire test was conducted with all equipment on two 
wood tables. One table was used for the magnet under 
test and the magnetometer probe while all other equipment 
sat on the second wood table. This was done to minimize 
magnetic interference due to ferrous parts in the test 
equipment. 
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Figure 7-3a 




Figure 7-3c 


Figure 7-3 Magnet Construction 
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II MAGNETIC ACQUISITION SYSTEM IMPLEMENTATION 

1.0 Introduction 

An automatic attitude acquisition system has been 
developed under contract NAS5-21649. A description 
of this system and its performance is given in report 
no. 90345, which is included as an appendix to this 
section of this report. 

The development of the hardware for this system was done 
under contract NAS5-21170, and is reported here. 

2.0 Description of hardware 

The hardware required for the magnetic acquisition system 
includes a 3 axis magnetometer, rate circuits, and 3 
magnets. The magnetometer developed under this program 
is similar to other magnetometers made by Schonstedt 
Instrument Company, but repackaged for the applications. 
On the following pages will be found: 

a) A preliminary specification for a three-axis 
magnetometer. 

b) A Schonstedt control drawing for the magnetometer 
built to this specification. 

c) Test data on the engineering model magnetometer. 

d) A photograph of the magnetometer probe. 
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1.0 SCOPE 

This specification establishes the interface requirements and 
performance of a set of three magnetometers for use in the SATS 
attitude control system. In addition, this document references 
documents which pertain to the quality of workmanship and parts 
to be used. 

2.0 APPLICABLE DOCUMENTS 

The following documents form a part of this specification to 
the extent specified herein. In the event of a conflict between 
documents referenced here and the detailed content of this speci- 
fication, the latter shall take precedence. 

2 . 1 NASA Publications 

NHB 5300.4 (3A) Requirements for Soldered 

May 1968 Electrical Connections 

PPL- 11 Revision A GSFC Preferred Parts List 

July 1971 

2 . 2 ITHACO Drawings / 

D41017 Frame, Electronics Univ. 

C88135 Printed Circuit Card, 

Control Drawing 

3.0 SPECIFICATIONS 

3.1 Electrical Characteristics 

The vendor shall supply three output voltages. Each 
voltage shall be proportional to one of three orthogonal 
components of the ambient magnetic field. A signal common 
shall also be supplied. The following specifications apply 
over the operating temperature range 0-50 °C. 

3.1.1 Operating Sources 

The magnetometer shall be capable of operating 
from externally supplied +10VDC and -10VDC ±1V. 

3.1.2 Sensitivity 

The three components of the magnetic field shall 
be related to the respective output voltages by 
the factor 4VDC/Oersted ±2%. 

3.1.3 Bi-polarity 

Reversal of the probe direction will cause the 
output voltage to reverse sign. 
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3.1.4 Minimum Saturation Levels 

The magnetometer shall not saturate when sub- 
jected to a magnetic intensity of 1 Oersted 
along the sensitive axis of the probe. 

3.1.5 Null 

In the absence of a magnetic field, the output 
voltage shall not exceed 20mVDC (at +25°G) . 

3.1.6 Linearity 

The magnetometer output shall be linear to 
within 0.5% of full scale. Full scale shall be 
defined as 0.80 Oersted. 

3.1.7 Ripple Voltage 

The maximum allowable ripple voltage at the out- 
put of the magnetometer is 20mV peak-to-peak. 

3.1.8 Temperature Stability / 

The change in null output voltage over the 
ambient range 0°C to 50 °C shall be no greater 
than 30mV. 

3.1.9 Noise 

After processing by a test circuit, the output noise 
voltage shall be no greater than 32mV peak-to-peak. 
The figure is a schematic of the test circuit. 

3.1.10 Crosstalk and Sensitivity to Crossfield 

With all three magnetometers energized, apply 
the magnetic field along the mechanical axis of 
One of the magnetometers. The resultant voltage 
in each of the other two aces shall be equivalent 
to not more than 1 degree in misalignment. The 
null voltage shall be subtracted before performing 

'- the calculations. 

\ \ 

3.1.11 Power Consumption 

Average power shall not exceed 750mV. 

3 . 2 Environmental Tests 

The launch and orbit conditions to which flight models of 
this design will be subjected shall be considered. 
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3.2.1 The satisfactory operation must be maintained 

over the temperature range from -10 °C to +60 °C. 
Some degradation in analog relationships is 
permissable at the temperature extremes, but 
all operations must function properly. 


3.2.2 The magnetometer shall be capable of meeting the 
requirements of Section 3.1 after being vibrated 
at the following levels with power applied. 

Sinusoidal Vibration 


Frequency 

Amplitude - "g" 0-to-Peak 

Range 

Thrust 

(Hz) 

Axis 

5-20 

0.5" D.A.D. 

20-30 

lOg 

30-85 

15g 

85-110 

37g 

110-2000 

15g 


Transverse Axis 

5-21 

0.5" D.A.D. 

21-30 

16g 

30-2000 

15g 


Random Vibration 


Direction 


Frequency 

Band 

(Hz) 


Spectral 

Density 

(g 2 /Hz) 


g-RMS 


Thrust Axis 20-2000 


0.2 


20 


Transverse 20-2000 0.2 

Axis 


20 


The duration of the test shall be 4 minutes in 
each direction — 12 minutes total. 


3. 3 Operability 

3.3.1 Maintainability 

The equipment shall be designed so that mainten- 
ance shall not be required during the specified 
operating and storage period. 
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3.3.2 Shelf Life 

The equipment shall be capable of meeting the 
requirements of Section 3.1 after a two-year 
storage period, under the following environments: 

Temperature 5 to 45°C 
Relative Humidity 50% or less 

3.3.3 Operating Life and Environment 

The magnetometer shall be designed to function 
in the free space environment for a period of 
three years. 

3.4 Circuit Packaging 

The magnetometer electronics shall be packaged subject to 
the constraints imposed by the drawings listed in Sec. 2.2. 

3. 5 Probe Packaging 

This section will be added later. 

/ 

3.6 Hand Soldering 

The provisions of NHB 5300.4 (3A) shall apply to all hand 
soldering operations. 

3. 7 Selection of Components 

Parts should be chosen from the GSFC Preferred Parts List 
(PPL-11, Revision A) . Parts which are not on the Preferred 
Parts List shall be approved by ITHACO. 

PREPARATION FOR DELIVERY 

The subcontractor shall package the magnetometer in such a way 
as to prevent damage or deterioration during shipment. Packaging 
shall be subject to review and approval by ITHACO. All shipping 
containers are to be plainly marked with ITHACO purchase order 
number, name, quantity, serial number and ITHACO specification 
number. 
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3.0 Magnetometer noise measurements 

One concern in the designs of this magnetic acquisition 
system was whether or not available magnetometers have 
sufficient sensitivity to be able to detect the very 
low fixed rates required. Three of the Schonstedt 
RAM5C-NB Magnetometers were purchased and evaluated, 
using circuitry similar to the differentiators that 
will be used in the final hardware. The results of these 
tests are given in the following report no. 90293. 
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MAGNETOMETER NOISE MEASUREMENTS 


PURPOSE OF THE STUDY 

This study is part of an attempt to detect apparent changes 
in the earth magnetic field from a satellite during initial 
attitude acquisition. A system designed to accomplish this 
(Figure 1) consists of the Schonstedt RAM-5C-NB Magnetometer 
followed by a differentiator circuit. The output of the differ- 
entiator would be fed to a threshold detector of an attitude 
control system. This report evaluates the noise type and levels 
present in the magnetometer-differentiator combination. 



<=,e OR- ^ KNv 


GENERAL OPERATION OF THE MAGNETOMETER-DIFFERENTIATOR 

The proposed magnetic aspect system (Figure 1) is composed 
of a probe and associated electronics (Magnetic Aspect Sensor) 
purchased from Schonstedt, and a differentiator circuit. The 
sensor provides a "d.c." output voltage proportional to the 
component of the magnetic field coaxial with the probe. The 
differentiator circuit not only differentiates, but provides 
attenuation over a range of frequencies outside of the range of 
vehicle rotation frequencies. Primarily, two circuits were used 
during the study, one with a corner frequency W c of 1 radian per 
second and the other with a corner frequency of 5 radians per 
second. The approximate equivalent circuits are shown in . 

Figure 2 . 
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The mathematical representation is 

Y(jW) = -jlOOW ? W is in radians per second 
(I+jtp2 
W c 

Figures 3a and 3b are the Bode plots corresponding to the circuits 
of Figures 2a and 2b. Zero db is defined equivalent to unity gain. 

TEST RESULTS 

Figure 4 is the actual circuit configuration of the differ- 
entiator with 1 radian/sec. corner frequency. Separate FET's 
were used. The circuit with a corner frequency of 5 radians/sec. 
employed the 2N3922 dual FET's. 


lo** 



subjected to each test. 

The outputs of the W c =l and W c =5 circuits were measured when 
these circuits were used in combination with magnetometer 3883. 
The two tests were not run simultaneously. The magnetometer was 
centered radially in a wooden cylinder which in turn was centered 
radially in an iron pipe. The pipe was then closed. A magnet 
external to the pipe supplied the magnetic field. The magnet 
produced a 2.3 volt magnetometer output. The output of the dif- 
ferentiators were recorded using a Brush recorder. 
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The peak to peak values obtained are 

8mV for the W c = 1 circuit, 

80mV for the W c = 5 circuit. 

Figures 5a and 5b are samples of the recordings from which the 
values were obtained. 

To determine the noise generated by the differentiator circuit, 
the output was measured with the input shorted. Both circuits were 
placed in a metal box and the outputs simultaneously measured using 
the Brush recorder. The peak to peak output was half a millivolt 
for the W c = 1 and 2 millivolts for the W c = 5 circuit. 

The magnetometers themselves do not have identical noise 
characteristics. Figure 6 is a Brush recording showing the noise 
output from the W c = 5 circuit fed by each of the three magnetometers. 
Each of the magnetometers was tested in the pipe with the magnet 
supplying the same field in each case. A considerable difference 
between the noise levels of the magnetometers is apparent. Dis- 
turbances, some deliberately created, in the laboratory verify that 
the magnetometers and the differentiator were functioning during 
the measurements. 

Two magnetometers, each in combination with a W c = 1 differ- 
entiator were tested simultaneously to show that spikes which were 
being observed on the oscioolgraph are not internally generated. 
Oscillographs obtained during the night hours also tend to confirm 
that the largest spikes are not internally generated. The night 
graphs are almost without the large spikes which are common in the 
data taken during normal working hours. 

Oscillographs were run for 1 hour with the W c = 1 differentiator 
and then for 1 hour with the W c = 5 differentiator. These recordings 
show that noise with a period on the order of minutes is negligible 
by comparison with the amplitude of the noise envelope. 

Root-mean-square (RMS) measurements were obtained to indicate 
the general level of the noise as a function of Wc. The measurements 
were obtained graphically from the Brush oscillographs. Three 
circuits were used in turn with magnetometer 3883. The magnetometer 
was enclosed in the pipe and the field supplied by a magnet external 
to the pipe. All three differentiators were of the same mathematical 
form, all had K = 100, and all had the same input capacitor and out- 
put resistor. With those restrictions, the noise levels may be 
considered to be a function of W c . The equivalent circuits for two 
of the differentiators are Figures 2a and 2b, the third (W c = 0.5) 
is shown by Figure 7 . 




Vertical Scale: lmV/small div. 
Horizontal Scale: 5mm/sec. 

Noise Output of Differentiator with W c = 1 rad/sec. in 
Combination with Magnetometer SN 3883. 

FIGURE 5A 



Vertical Scale: 5mV/small div. 
Horizontal Scale: 5mm/sec. 

Noise Output of Differentiator with W c = 5 rad/sec. in 
Combination with Magnetometer SN 3883. 

FIGURE 5B 



3882 3881 3883 

-N_ 

Comparison of Noise Levels from Three Magnetometers. 
Output of Differentiator with W c = 5 rad/sec. 

FIGURE 6 
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of the graphical analysis are given in Table 


1. 


W c 

radian/sec. 


Magnetometer 3883 Output Differentiator Output Noise 
Volts d.c. Millivolts RMS 


0.5 1.6 0.6 

1.0 2.3 1.3 

5.0 2.3 8.1 

TABLE 1 

To obtain the RMS voltage graphically, the instantaneous value 
of the noise voltage is measured at 1 second intervals. The square 
root of the mean of the square of these measurements is the RMS 
voltage. 

MAGNETOMETER OUTPUT CHARACTERISTICS 

Tests show that the d.c. output impedance of the magnetometer 
is 3.6K ohms. The 20K ohm resistance shown in Figure 2b actually 
includes the 3.6K ohm output impedance. 

At a frequency of W c = 0.5 radian/sec., the impedance of the 
10y farad differentiator input capacitor is 200K ohms, so there is 
neglibible loading from d.c. to W = 0.5 radian/sec. 

Discussions with the vendor indicate that no electrical isolation 
problems exist in connection with the output of the magnetometers. 

CONCLUSIONS 

To get some idea of the feasibility of using the particular 
magnetometer-differentiator combination studied, some assumptions 
of physical conditions will be made. 

For purposes of this study, a satellite will be assumed in a 
500 mile high circular orbit. The RAM-5C magnetometer has a 
sensitivity of 4 volts per oersted. Assuming a maximum sea level 
field strength of 580 milliguass in the orbit plane, and an inverse 
cube variation of field strength with orbit radius? a voltage of^ 

1.6 volts maximum would be present at 500 miles. 

The Bode plots. Figure 3, give the voltage gain of the corre- 
sponding differentiators as a function of the frequency of a sinu- 
soidal input. For simplicity, assume the satellite to be rotating 
about an axis which is fixed relative to the earth. Assume we are 
interested in rotation frequencies from 0.5 radians/sec. down to 
the orbit rate, and assume the orbit frequency is .001 radian/sec. 
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According to the Bode plots, a magnetometer signal at the orbit 
frequency would be attenuated by a minimum of 20db. For the 500 
mile orbit, the peak differentiator output would be a maximum of 
±160mV. If the angle between the magnetometer probe and the orbit 
plane were such as to reduce the field strength by a factor of 0.1, 
then the signal would be attenuated another 20db and the peak output 
of the differentiator would be ±16mV. 

Suppose threshold detectors to be at the output of the dif- 
ferentiator. Based on the above discussion, the threshold levels 
might be set between ±16mV and ±160mV. Recall that the peak to 
peak noise values recorded on the oscillograph were 

8mV for the W c = 1 circuit, 

80mV for the W c = 5 circuit. 

The oscillographs show that the noise level measured at the 
output of the differentiator is dependent on the d.c. output voltage 
of the magnetometer. The reason is not known. 

The study does not conclusively define the limitations of this 
particular magnetometer-differentiator combination. However, based 
on the measurements , the prospects for this combination appear 
promising. 
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MAGNETOMETER NOISE MEASUREMENTS 
Addendum to Report No. 90293, Aug. 27, 1971 

Report No. 90293 gives the results of noise measurements 
for the Schonstedt RAM-5C-NB Magnetometer in combination with 
a differentiator circuit. As mentioned in that report, the 
noise level measured at the output of the differentiator is 
dependent on the d.c. output voltage of the magnetometer. 

Oscillographs of the noise at the output of the magneto- 
meter itself shows a dependence on the d.c. output voltage of 
the magnetometer. To obtain this more recent data, magneto- 
meter 3883 was enclosed in the iron pipe and a constant magnetic 
field supplied by a magnet external to the pipe. The output of 
the magnetometer was coupled to the preamp of the Brush Recorder 
through a lOyf capacitor. Samples of these recordings are pre- 
sented in the Figure. 

A close look at the oscillographs show the noise level to 
be higher for a magnetometer output voltage of -2.4V d.c. than 
for an output voltage of +2.4V d.c. Since all of the data in 
Report No. 90293 was obtained for positive output voltages, the 
following question arises. If data were taken for negative and 
positive instead of just positive output voltages, might magne- 
tometer 3881 be found to have greater p-p noise than 3883, 
instead of the other way around. To resolve this doubt, 3881 
was tested in combination with the Wc = 5 differentiator; the 
magnetometer was enclosed in the pipe and the field supplied by 
a magnet external to the pipe. 

The following table shows the p-p differentiator output 
noise as a function of magnetometer output voltage. The data 
is only approximate due to the difficulty of defining p-p noise. 
However, this data clearly shows that the noise levels are well 
below the 80mV p-p measured with 3883. 

Differentiator Output Noise vs Magnetometer 3881 DC Output 

Magnetometer Output W c =5 Differentiator Noise Output 

Volts DC Millivolts p-p 

-2.4 
- 2.0 
- 1.0 
0.0 
+ 1.0 
+ 2.0 

+2.4 35 


25 

25 

14 

9 

20 

19 
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CONCLUSIONS 

Noise levels vary from magnetometer to magnetometer. 

Noise levels also vary with field strength. Finally, and most 
important, if rates on the order of orbit rates are to be de- 
tected, the results of the tests with the W c = 5 differentiator 
are inconclusive. Simulation will be needed to provide the 
answer. However, if a relay were installed to convert the 
W c = 5 differentiator to a W c = 1 differentiator during acqui- 
sition, detection at orbit rates should be feasible. 


Jay Levy 
ITHACO, Inc. 

735 W. Clinton Street 
Ithaca, New York 14850 
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ABSTRACT 

A Magnetic Control System for Attitude Acquisition 


This report describes a spacecraft magnetic attitude 
acquisition system that is capable of automatically despinning 
a satellite from arbitrarily high rates around any axis, and 
provides terminal orientation that makes capture by conventional 
fine control attitude control systems routine. The system con- 
sists of a 3-axis magnetometer, a set of 3 orthogonal magnets, 
and appropriate cpntrol logic. No earth sensor is required. 

Acquisition is treated in two phases. During the despin phase 
we are concerned with removing the tumbling motion of the satellite 
In this phase the performance of the system is unaffected by the 
presence or absence of a momentum wheel. Phase 2, that of 
orienting the spacecraft to the desired attitude, requires that a 
momentum bias wheel be present. In the terminal orientation, the 
axis of the momentum wheel is substantially perpendicular to the 
orbit plane (roll and yaw errors near zero) and the pitch rate of 
the satellite is at twice orbit rate. 

This report describes the analysis and simulation that has 
been done in evaluating the performance of this system. A well- 
configured system will result in despin times of the order of 5 
orbits per RPM for spacecraft in low earth orbits. Following 
despin, terminal orientation is achieved after another one to 
three orbits, depending on the capture range of the associated 
fine control system. 

While this report does not describe the physical hardware, 
the system can be implemented inexpensively with weight less than 
5 lbs. and power of about 3 watts. 
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INTRODUCTION 
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In this report we will -be concerned with the description and 
evaluation of a magnetic acquisition/attitude control system. The 
work described here has been performed during the last several 
months, principally by the author, with the assistance and guidance 
of Messrs. David Sonnabend, and Robert Fowler. This work was done 
under contract #NAS5-21649, at the behest of the Systems Analysis 
Branch, GSFC, Greenbelt, Maryland. The original concept of this 
system was suggested by Seymor Kant, Head, Systems Analysis Branch, 
GSFC. 


The attitude control system described here is designed to 
"despin" a satellite and then orient it to any preferred Earth 
referenced attitude. This is done thru interaction with the 
Earth's magnetic field. An interesting feature of the system is 
that it requires no Earth sensor. 

In this report we will discuss the control law for despinning 
the satellite and the mechanism thru which some preferred attitude 
is attained. Then we will discuss stability considerations as they 
apply to this system. Next, system performance is evaluated in the 
light of both analytical results and simulation results. We will 
be concerned with both despin and attitude acquisition performance. 
Finally, some recommendations for further study are included in 
case any parties might wish to continue the work discussed here. 

Appendix B displays a listing of the FORTRAN program used in 
the digital simulation. Appendix C reviews some of the dynamic 
stability theory for dual spin satellites, for convenient reference, 
but does so in a different fashion than is usually employed. The 
results obtained are, of course, the same as usual. In Appendix D 
we analyze the equilibrium points of our system; the stability of 
these points is briefly treated in Appendix E, using some of the 
results of Appendix C. 

The author wishes to make one further remark before continuing. 
The (differential) equations describing the behavior of the system 
discussed herein are nonlinear and otherwise intractable. This 
being the case, it is not possible to obtain explicit analytical 
solutions and recourse to simulation must be made. Much of the re- 
sults presented here were obtained via a digital simulation model. 

A program, previously written by the author and noted in the refer- 
ences (1) , was used to this end. A listing of this program, and an 
index of computer runs, may be found in the appendices. A brief 
description of the simulated satellite is contained in the following 
section . 


1 


f 
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SIMULATED SATELLITE 

Most of the simulation runs dealt with a satellite with the 
following characteristics: (1) moments of inertia of 47.5Kg-m 2 

(35 slug- ft 2 ) , 67.9Kg-m2 (50 slug-ft 2 ), and 33.9Kg-m 2 (25 slug-ft^) 
about, respectively, its yaw, roll and pitch axes; (2) a pitch 
momentum bias of 0.944 Nt-m-sec (0.7 lbf-ft-sec) ; and (3) switch- 
able control magnets of 10^ pole-cm strength. There are three such 
magnets, one each along the yaw, roll and pitch axes. We again 
note that these are principal axes for the satellite. This satel-r 
lite is in a circular low altitude orbit with a period of approxi- 
mately 100 min. We note however, that these parameters of orbit 
and vehicle are of no special significance, but were chosen merely 
to implement the digital simulation. 

DESCRIPTION OF PROBLEM 

In general, a satellite injected into orbit will be observed 
to be tumbling about with a certain residual angular velocity. This 
is so in spite of efforts made to release it "gently". Various de- 
vices have been utilized to dissipate this residual motion. In our 
case, the satellite will be launched by a spin stabilized rocket and 
hence will be inserted into orbit with an angular velocity of roughly 
100 RPM about some (known) axis. A simple mechanism (a yo-yo) will 
then reduce that spin to the order of 1-10 RPM. However, this rate 
is, in general, still so high as to prevent the basic attitude con- 
trol system from locking onto some target (the Earth, sun, stars, 
etc.). A further spin rate reduction is required, preferably to say 
0.1 RPM or less. To this end, such active devices as rate gyros, 
accelerometers, and magnetometers coupled to gas jet systems, and 
passive devices such as eddy current and magnetic hysteresis rods 
have been employed. 

Reverting to general terms, then, there are two distinct 
phases of the attitude acquisition procedure. The first is despin, 
of which we have just spoken, and the second is an attitude orien- 
tation phase which facilitates capture by the basic control system. 
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EQUATIONS OF MOTION 

The equations describing the rotary motion of a satellite 
(or any other rigid body) are most conveniently written for a set 
of body fixed principal axes. They are the well known Euler's 
Equations of Motion, with terms added to account for momentum 
wheels. To wit, 


t = d^ PACE (H) = d^ 0DY (H) + a) x H ( 1 ) 

dt dt 

where H = I • or + h ( 1 ') 

or, in scalar form 

= h^ + + t02 h 3 — ^2^2 + ^ 2 W 3^3 — ^2^ (l a ) 

T 2 = ^2 + l 2 w 2 ■*" W3hj_ - o)]_h.3 + (133(11 ~ I3) (lb) 

x 3 = h 3 + I 30)3 + w^h2 ~ ^2^1 w l w 2 ^2 ~ i]_) (lc) 


The reader is cautioned that great care should be exercised if 
these axes are not right handed. As mentioned above, these axes 
are principal axes passing thru the center of mass. Here,_the Ii 
are the moments of inertia about the ith coordinate axis, t is the 
applied torque, and h is the momentum contribution of the wheels. 
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CONTROL THEORY 

As noted above, the attitude acquisition system performs two 
principal functions. The first is despinning the tumbling satel- 
lite, and the second consists of providing some preferred orienta- 
tion. Consider the former. 

The rotational kinetic energy" of the satellite, not including 
that of any reaction wheels, is 

2 

T = 1 Z IiUjL (3) 

2 i 

We would like to diminish this quantity. 

d (T) = T = = x * to (4) 

dt i 

where t is the external torque acting on the satellite, and 
t • a) is the _rate at which work is done on the satellite. 

For a dipole M in a magnetic field B, the torque exerted by the 
field on the dipole is 


T = M x B 

Combining (4) & (5) 


(5) 


T - M x B • to = B x to 
Now consider the quantity 

df ACE (B) = df DY (B) 
dt dt 


d_ 

dt 


M 
BODY 


(B) 


+ oo X B 


( 6 ) 


(7) 


SPACE - - 

Suppose d_ " J (B)=0. That is, we say that this term is negligible 
dt 

in comparison with the other two terms. Then (7) becomes 

d^ 0DY (B) = B x w (8) 

dt 


Substituting in (6) we obtain 

T = B • M (9) 

2 , The desired control scheme is now clear. We must simply measure 
B along any axis in the satellite, and then change_the polarity of a 
magnet lying along that axis, keeping the sense of M opposite that of 
«£» t t • ( t ( 

B. This insures T < 0 and decreasing T. The limitations of such a 
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scheme are also clear. Equation (9) is valid only when (8) is 
valid, that is, if d SPACE (B) is truly negligible in comparison 

_dt 

with a) x B. Now, B varies in inertia space_as sin2w 0 t, where 
oj 0 is the orbital frequency. Thus, d SPACE (B) is equivalent to 

-dip 

a body rate of 2to 0 / and as long as w is, say, an order of magni- 
tude greater than 2w 0 (about 0.002 rad/sec for low altitude 
orbits), we may safely neglect d SPACE (B) . 

dt 

For reasons which will become clear later, the chosen con- 
trol configuration involves three magnets. These lie along the 
principal axes of the vehicle and are controlled by magnetometers 
(which measure the field component) along these same axes. The 
controller operates in a flip-flop manner, switching the magnets' 
polarity to keep the M-j_ opposite in sense to fei , as measured by 
magnetometers on board the satellite. 



Report No. 90345 
Page 6 


TERMINAL ORIENTATION 

We now consider the second phase of attitude acquisition, 
orientation. We wish to attain some particular attitude. As long 
as this preferred orientation is with respect to the Earth, any 
desired attitude may be attained. In this section we will confine 
our discussion to a satellite having only a pitch momentum bias. 
Other momentum wheels, if any, are not yet activated at this point. 
The arguments advanced below are not as rigorous as might be de- 
sired but it is felt that their physical appeal more than compen- 
sates for that. 

Consider the satellite to be orbiting the Earth at angular 
rate w 0 , and assume also that it has been completely despun. Then 
it must be that the satellite interprets this situation as though 
it were spinning backwards at angular velocity -2u) 0 e n , where e n is 
a unit vector normal to the plane of orbit (and whose sense is 
determined by the direction of orbit) . This is because the satel- 
lite senses the Earth's magnetic field passing by it twice per 
orbit. Note that the trajectory yaw and roll axis components of 
the Earth's field vary as sin 20, where 0 is the angle into orbit 
from some reference point. Since the system "senses" an angular 
velocity whose direction is -e n , it responds with a corrective 
torque along the direction +i n » The net result of such a torque 
must be to gradually align the system momentum bias (which properly 
lies along the vehicle pitch axis) with the orbit normal. In this 
way the proper yaw- roll attitude is attained. It is the fact that 
the satellite is not spinning at 2wo that results in the net torque 
along e n which in turn aligns the momentum bias with the orbit 
normal. But what if it had been spinning at 2w 0 originally, or in 
some way reached this state. It is apparent that this could occur 
in only two ways. That is, the system momentum must be either 
parallel or antiparallel to e n . If this is not so, then a component 
of this momentum lies in the orbit plane and a body fixed torque, 
also lying in the orbit plane and of the correct magnitude would be 
required to force the body to rotate about e n at 2coo- Where is this 
torque to come from? The author does not believe in its existence 
and hence concludes_that the system momentum bias must lie parallel 
or antiparallel to e n if it is to rotate at 2io 0 . Both of these 
situations correspond to equilibrium points. We simply note here 
that the antiparallel situation is not a stable condition. 

Our conclusions ,are then that the satellite will eventually 
align its momentum bias axis (in our case, the pitch axis) with the 
orbit normal and thus provide yaw- roll attitude acquisition. If 
the reader is not completely convinced of this by the rather un- 
rigorous arguments advanced here, we reassure him by noting that 
considerable simulation experience is completely in accord with the 
above analysis. A further discussion concerning the amount of 
momentum bias required in the above scheme and the effects of more 
than one momentum wheel, etc., are reserved for later sections. 

Refer to the section titled "System Performance" and see Appendix C. 
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STABILITY CONSIDERATIONS 

The desired action of the control law just developed is, 
of course, to reduce the angular velocity of the satellite to 
zero. Such a condition is not, properly, an equilibrium point 
for our system, since no provision has been made, as of yet, for 
a null state for the magnets. They are always on, having one 
polarity or the other, and the satellite always has some angular 
acceleration, hence the system possesses no equilibrium points. 

It would be easy however, to insert deadzones in the magnet 
controllers. Provided these deadzones are large enough to_prevent 
magnet activation due to the relatively small term dSPACE (§) in 

dt 

(7) (which the satellite experiences due to its orbital_motion and 
the rotation of the Earth) , then, the condition to = 0, B = any, is 
an equilibrium point. If this sort of provision is not made, then 
it would seem that the satellite should behave in some sort of 
quasi-oscillatory manner, at small angular rates. That this con- 
dition would, once attained, persist, seems quite plausible from 
the analysis above. It is not, however, a simple matter to talk 
about equilibrium points and stability. The equations describing 
this system are quite intractable to ordinary stability analysis. 
For one thing, the control law is binary. Even if we replace it 
with a linear law (the first term of its describing function) the 
resulting system equations consist of six first order nonlinear 
differential equations; three involve quadratic combinations, the 
other three cubic combinations of the state variables (wi, o> 2 , wo, 
B lf B 2 , B 3 ) . 

We should note here that if the system has either a binary 
control law with deadzones or a linear control law (that is, M = 

-KB) , then there are an infinity of proper equilibrium points. 

These points belong to one of four distinct classes. However, in 
order to spare the reader, a discussion of these points and their 
stability is relegated to Appendix D. Suffice it to say that 
some of these points (other than w = 0) can be shown, by analysis, 
to be neutrally stable. Further, simulation results show that in 
some domain some of the points are stable. 

There are other problems too. We are particularly interested 
in investigating the existence of any so called "psuedo equilibrium 
points" (PEP's for brevity). Since dSPACE ( b) 0, and there are 

J 4- 

various disturbance torques acting on our satellite, and the 
system equations are nonlinear, there may be points which are not 
properly equilibrium (stationary) points but near which the system 
might "hang up". Here we are talking principally about limit cycle 
behavior, but must also be concerned with various (as yet unknown) 
forcing functions. 

In view of these difficulties, it was decided that digital 
simulation was the most suitable tool for evaluating this system. 

A discussion of the results of this simulation work is reserved 
for the next section and Appendix E. 
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Although the actual stability analysis is put off until 
appendices C and E, the conclusions and implications of those 
analyses are reviewed here. First, we note that, for a vehicle 
configured as described in the section "Simulated Satellite" 

(that is to say, the proposed vehicle), the terminal condition 
(or behavior) corresponds to a stable equilibrium point. This 
is, of course, most desirable. We again note that the "upside 
down" orientation (that is, with the pitch axis flipped over) 
corresponds to an unstable equilibrium state; this too is 
soothing. Finally we note that, for our vehicle, a spin rate 
of up to 27.8 x 10 - 3 radians/sec. about the pitch axis also may 
be a stable equilibrium state. 

This stability conclusion neglects the fact that the direc- 
tion of B in space changes as the satellite orbits the Earth. 

This turns out to be a saving grace. Nonetheless, for strong 
enough magnets, the satellite may "track" the field. During such 
a period it is possible that the despin rate may be reduced by as 
much as a factor of fifty. For this reason it might be advisable 
to insert a mode switch for disabling the despin system for an 
eighth orbit or so if_it seems as though this has occured. At 
the end of that time w and B would no longer be parallel and hence 
the despin system will again perform properly. 
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SYSTEM PERFORMANCE 

Having convinced the reader (hopefully) that our system works, 
we now attempt to answer the question - how well? As before, we 
shall first be concerned with the despin phase. 

As the system kinetic energy j wl*w tends towards zero (this 

was assured in the section titled "Control Theory"), so must to and 
the system momentum j. % (this, of course, excludes the wheels' 
momenta). j n fact, the despin process may be looked on as a re- 
duction of the momentum term to zero. In the very best case 
this requires a time At op t given by 

Atopt = 1^1 = [Pol (10) 

T max MB 

Here P Q = is the initial momentum of the body other than 

that due to the wheels. It is useful to consider an "efficiency 
factor" a, 

a = Atopt (11) 

Atactual 

where a of course lies between zero and one. a is a function of 
various parameters, such as the inertia configuration of 

the satellite, and the altitude and inclination of_orbit. That a 
is less than one is a reflection of the fact that M is not always 
perpendicular to B and that t does not always lie opposite 

For a given system configuration (I & M) and a given orbit 
and initial conditions (determining to, B) we should be able to de- 
termine a and hence the time for despin. This was the goal of much 
of our early simulation, in which various orbits, initial conditions 
and system configurations were evaluated. This simulation shows 
that a generally lies in the range 0.55-0.80. For the worst case, 
a near equatorial orbit with initial spin along the field vector, 
a may be as poor as 0.15-0.20. 

Another way to look at the performance of this system is from 
an energy point of view. Consider, for simplicity rotation about 
a single axis. 

T = jlu 2 (12) 

T = 1 0000 = WT (13) 

• * Iw = t = aMB (0<a<l) (14) 

The last of these equations says that w, and the angular momentum 
Iu), are always decreasing. The constancy of this rate of decrease 
depends on the constancy of a. We are of course assuming that t 
has the proper sense. 
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The attitude acquisition part of the problem is more difficult 
to analyze. If, for example h approaches 0 we will not _ acquire. 

But also if h approaches infinity we will acquire, but it will t k 
infinite time. Also, if M approaches 0 or infinity we will not^ 
attain the proper attitude. In view of these consideration 
seem that there might be some ratio of momentum bias to magnet 
torque which would provide optimum performance. By this we mean^ 
both minimum time to acquire and minimum roll-yaw error after ac 
quisition. Consider the parameter 


T S ys = _h_ 
MB 


(15) 


which has the dimensions of time. It is rather difficult to say 
anything about this problem analytically, but the following is 
observed from simulation runs. 

We first note that there is a certain minimum bias required for 
terminal attitude acquisition. It is shown in Appendix C, eq. (CIO) , 
that 


(CIO) 

2 wo (for t he terminal condition) .; 

(16) 


h 0 > max (Ij. , 12 ) 

for us, ho = h + I 3 S and (2 
Thus (CIO) becomes 

h > 2 (do { max ( 11 , 12 ) ” I 3 ) 

For good performance this amount should be exceeded by a factor of 
four or more. Little is gained by increasing h beyond this point. 

If we now choose a particular orbit, the time and quality of ac 
quisition (after despin has been completed) are functions of the 
magnet strength M. Acquisition time is also, of course, very much a 
function of attitude at the conclusion of despin. A (very) quanta 
tive picture of this situation is presented m Figure 1. As pointed 
out previously, for very weak or very strong magnets we do not attain 
the desired attitude. The author cannot caution the reader too 
strongly against taking this curve too literally. _ Much more work 
would be needed to get a more exact picture of this phenomena. There 
are too many variables and too few points on the curve of Figure - 
to let the author rest comfortably. As for the quality of roll yaw 
acquisition, it is interesting to note that a curve qualitatively 
similar to Figure 1 may be used to describe this as well. Figure 2 
depicts this situation. Note that the minimum point is shifted to 
about 2000 secs. (This appears to be significant) and the curve is 
somewhat flatter. The best, that we can do appears to be about 2 . 
Again, the author cannot over-emphasize the tentative nature of these 
findings. Considerable work remains to be done in this area. 

For a well proportioned system, acquisition times are of the 
order of one orbit. This being the case, it is expected that ^ 
optimum performance will generally be defined in terms of quality 
of roll-yaw acquisition rather than in terms of minimum time to 

acquire. 
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RECOMMENDATIONS FOR FURTHER WORK 

The author feels that the work done to date definitely demon- 
strates the feasibility and desirability of this acquisition con- 
trol system. Nevertheless, work remains to be done. 

There is the matter of system stability. While it seems 
apparent that there are no stability problems which might interfere 
with the proper operation of the control system, nonetheless it 
would be good if we obtained some further analytical results. In 
particular, an analysis including the effects of the change in B 
due to orbital motion would be interesting. 

As noted, we have not included the effects of disturbance 
torques on the performance of this system. Their effects on both 
limit cycle behavior and terminal attitude orientation should be 
investigated. Their neglect has been justified up to this point 
by the fact that the magnet torques are at least an order of mag- 
nitude greater than the disturbances. 

Another area which needs refinement is the magnetic field model 
used in our simulation. So far, only the simple tilted dipole has 
been used. The effect of a more accurate field model on terminal 
attitude orientation needs to be investigated. 

More investigation into the effects of deadbands, hysteresis 
loops, and inhibit controls placed in the magnet control loops would 
be desirable. There is also the matter of crosstalk between magnets 
and magnetometers. We must be sure that no unstable loop is set up 
in which switching magnets trigger magnetometers which in turn cause 
the magnets to switch again. 

One other point which will bear further investigation is the 
potential of using this acquisition system as a complete attitude 
pointing system. Used with a pitch reaction wheel scanner and simple 
pitch control loop we would have complete three axis control. Some 
studies of pointing accuracies versus various system parameters would 
be useful in investigating this possibility. 
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DIGITAL SIMULATION PROGRAM LISTING 



Report No. 90345 


This appendix consists of a listing of the digital computer 
program used to simulate control system/satellite behavior. The 
program is quite similar to the one discussed in reference 1, and 
the reader is referred there for a more detailed discussion, flow- 
charts, etc. Its operation may be briefly summarized as follows. 

The main program performs the bookkeeping functions of input/ 
output and maintains overall control. What we are doing here is, 
essentially, integrating equation XI) , specialized for the vehicle 
described under "Simulated Satellite". This is done through sub- 
routines WDOT and RK4 . An additional problem we must deal with 
here is thht the applied torques are functions of the vehicle's 
attitude. Namely, they are functions of the Earth's magnetic field 
components as measured in the vehicle. For this reason we must 
keep track of the vehicle's attitude. This is implemented through, 
the use of Euler Homogeneous parameters in subroutines . EHPQT , NORM, 
MATl, and MULVEC. The Earth's magnetic field is first obtained 
in trajectory axis components (subroutine BFLDS1) , and then trans- 
ferred into body axis components, making use of the direction 
cosine matrix (subroutine EHPA) . Subroutine TORQUE then computes 
the torques arising from the interaction between this field and 
the on board magnets. Please note that the subroutine displayed 
contains neither deadzones nor hysteresis loops, although some 
subroutines incorporating them have been used at times. / 


B-l 



FORM NO '.in:* PkSMll'O 


-JOB 

c ****CE MfRAL PURPOSE KINEMATIC INTEGRATOR PACKAGE USES El/LER DAV20001 

C ~*#*^D!WGErfEaUS PARAMETERS IN ORDER TO SOLVE PROBLEMS INVOLVING DAV20002 

c ****booy kate specified -rotation rates® ******* ************ da v 2 0003 

c ****ftUTHOR: ALBERT C* STICKLER DAV20004 

q *#**: *#**><<##*#>:<**:**********:*#************************************* ****DAV2C00 5 

C~ « w H AIN PROGRAM SKELETON FILL IN DETAILS," PARAMETERS , 'ETC* *>T****DAV2 0006 

C **** ASSEMBLAGE INCLUDES SUBROUTINES NORM, MAT l » ~RK4 t EHPA , _DAV 200.07 

C ****EH'PgtV SANOLEt AND INTERP , AMONG OTHERS ************************OAV20008 

C DAV20009 


1 

2 

REAL 11, 12, 13, QT(4,A), A(3,3)»AXIS(3) 
REAL*** INCL lN»“BT(3 ) 


DAV20010 
DAV200 11 


3 

4 

COMMON m/P C4) , W(3) ,QTfVTf WO 

C OMMON /WDO TS/ TC3) , C OWN (10), Cl, C2, C3, C4, C 5, C6 


DAV20012 

DAV20013 


5 

6 

REAL T3 ME/0 . 0/ , THE TA/0„ 0/ , WE /0 o 41 784E-02/ 
INTEGERS I/-1/, THEAW-l/ 


DAV200 14 
DAV20015 


7 

READ ( 5 , 90 ) 11 , X2 ,X3 , BIAS , FAC TOR 


DAV20016 


8 

READ (5, 30) PRINT, BSTEP 


DAV20017 


9 

READ (5 ,90)W 


DAV20018 


10 

■ReAPC5,90)I>T,rRUN 


DAV20019 


LI 

■READ t5, 90>1NC LIN, EARTHR , ALf 


DAV20020 


12 

■PEP IQP=0« 3 151243* SQRTCO • 1 *ALT ) **3 


DAV2002I 


13 

WRITEtE, 110) Tl, 12, T3, 3 IAS, FACTOR 


DAV20022 


14 

WRITEU,ltU -PRINT, BSTgp 


DAV20023 


15 

WRITE (4,112) DT, I NCL I N, EAR THR , ALT, PERIOD 


DAV20024 


16 . 

C T er CI2-T3i /II 


DAV20025 


17 

c2=-d3»ia )r& 


DAV20026 


13 

C3s(Xl -T2 ) /I.3 


DAV20027 


19 

e4= 1/.T1 


DAV 200 28 


20 

esn/x2 


DAV20029 


21 

C6=i/I3 


DAV20030 

►8 w 


c REAU AM) WRITE initial attitude 


DAV20031 

(D (D 
iQ TJ 

22 

READ(S,40)P 


DAV20032 

(D 0 

H 

23 

write ( 6 , 9 n p 


DAV20033 

tt) rt 

24 

~91 roAMAT r/ / , 1 INITIAL EHPS WERE 1 , 4F10a 4) 


DAV20034 

NJ 2 

25 

CONST = 0*00 1* FACTO R 


DAV20035 

O 

• 

26 

C0MMK 1 ) =#3A^C4 


DAV20036 


27 

C0MM(2 )=’8IAS*C5 


DAV20037 

"o 

28 

X{>STCT=BSTEP/360. 0*P£Jt lOD/DT 


DAV20038 

U) 

29 

XPjH*t7 = PRlNT 


DAV2003S 5 - ’ 

Vi 

30 

PTHETA = 360. 0/PERXQD*0T 


DAV20040 


31 

C At L 3FLDS 1 ft'NC LXN , ALT ,3T ) 


DAV20041 


32 

CALL BFLDS2IEARTHR , THETA, BT) 


DAV20042 


33 

WO = 6 2832/PERIOD • 


DAV20043 



■OrWAFsDT^P , 5 

, IZlTolt 

IF(jX/IPRIMT«'IPRINT .HE. II GO TO 75 

ERHOR=PMAC-1.6 ■ — UAV2L047 

ENERBWK’Wi 1 >«W( I H-I2»W(2 )*HI2H.I3*WI7)«III^I 

IHfAO=IHEA»+l nsi,,nr>.'n 

i r un rAP/50 «50 .eo. iHEad i write i«, ion 

^7^=1 t^Sr TI " E *^ r, ®5 '^RRO^VERERGV f A( 1 ,3) , A(2 , 3 1 HlWfis ! T 

CAU CHRA <f>,A) — o? JJ? 

theta=th£t a+ ptheta DAvloo« 

rc/cOcoiw^f-? S ^"^ ’ XICALL ^FLDS2 ( EARTHR + WE*T 1 ME tTHETA , BT I 0AV20056’ 

JFIg^EPGT .L T. 0.01) CALL BFLPS3 i THETA, BT > 

CALL MULV£C(A,BT,BB> — ®V>6oir 

CALL TORQUE I T .BB , CONS n DAV20059 

CALrRKA " DAV20060 

T F ( TJM E «€ T* TRU N ) STOP DAV20061 

GO TO 50. DAV20062 

90 FORMAT (5F1Q«3) DAV20063 

100 FORMAT (IH ,F6,6, 2 X ,3 ( £J0 o 3 , 1 X ) , 3 ( F6. 3 , 1 X ) , 2E 1 5 : 5 , 2F7o 3 ) DAV2006A" 

101 FORMAT ( IHl , 1 TIME ONEGA i 0MEGA2 0ME&A3 B1 BDAV20065 

62 B3 ERROR ENERGY +ROLL NES, W r > DAV20066 

1 100 F0R MAT( ///// , 50X, 1 SYSTEM PARAMETERS ARE ',//, 60X , 'INERTIA ( 1 >- • , DAV20067 

L FI 0. 2 » / / ,60 X , • INERTIA (2 >-• , FlO«2 ,7 / t60X, * I NERTI A( 3 ) = ' , F10 P 2 , ““ DAV20068 

?L //f 60X, 'MOMENTUM fllAS^' , FI 0,3, //,60X , » MAGNET STRENGTH FACTOR=* ,DAV20069 

3 F 1 0 . 2 ) — D A V 700 70~ 

1 1 10 FORM AT (///// 1 50 X I NTERVALS ARE' ,// ,60X, » PR INT EVERY* ,F10oO, DAV20071 

1 ' TIME STEPS • i f 7 ,2 5Xi 'CALCULATE MAGNETIC FIELD EVERY 1 , Fl'Oo'2 r~DAV2C072 

2 _ ' DEGREES OF ORBIT') DAV20073 

1 12 OFO RM AT (77///, 1 0 X , ' JnTeGRATI ON STEP SIZE=' , FS.2,' SEC ONDS , -- RUNGPAV2 0 077 

LE-KUT TA FOURTH ORDER INTEGRATION WITH EULER HONOGEN EOU $ TAR AMETERSDAV 20075 

2 ',//// ,50X,'0R6IT SPECIFICATIONS' , //,25X, 'INCLINATION ANG* = ' , F10, 2PAV20076 
3 . .. * 1 DEGREES' ,//,25X, '.ROTATION OF EARTH MINUS LONGITUDE OF THE. DAV20077 

4 MODES (INITIAL )= ' »Fl0o2, • DEGREES' , // ,50X, • ORBIT ALTiTUpE < CENDAV20078' 

STRAL DIST ANCE )=' ,F10»0,« KM* //,25X, 'PERIOD FOR THIS ORBIT= • , DAV20079 

6 FlOaO, • S£CS<* 'iff / ,25’X , ' ALL UNITS ARE IN M-K-S S Y S T E M~UNL E S S DAV20080 

7 ,,_°THEaHISE NOTED') DAV20081 

ENJ) DAV20082 


9|0 34 


60 ‘ SUBROUTINE EHPA(P»A) " " DAV20083 

61 REAL*4 P(4), A { 3 » 3 ) DAV20084 

C ****THIS ROUTINE RETURNS THE ROTATION MATRIX A CORRESPONDING DAV20085 

C ****TO THE El U L ER HOMOGENEOUS PARAMETERS P SUB I a A* V ( OLD ) = V( NEW ) DAV 2 00 86 


62 

63 


P 1 2 = P ( 1 )*P ( 1 ) 
• P22=P(2J*P<2> 


DAV20087 

DAV20088 


64 

65 


P32=P( 3) *P ( 3) 
P42 =P ( 4 ) *P ( 4 ) 


DAV20089 

DAV20090 

! 

66 

67 


A ( 1 , 1)=P42+P12-P22-P32 
A ( 1 , 2 ) =2o * ( p U. ) *P •: 2 ) + P( 3 ) *P( 4 ) ) 


DAV2009 1 
DAV20092 

! 

68 


A(l,3)=2o*(P(l)*P(3)-P(2)*P(4)) 


DAV20093 


69 


A(2,l)=2 3 *<P<l)*Pf2)-P(3>*P{4)> 


DAV20094 


70 


A (2, 2) =P42+P22-P32-P12 


DAV20095 

{ 

71 


A( 2,3 ) =2c * ( P(2 )*P( 3) +P(1 >*P( 4) ) 


DAV20096 


72 


A< 3* i)»2o *( P( 1)*P(3)+P( 2)*P(4) ) 


DAV20097 


73 


A( 3,2)=2 0 *(P(2)*P(3)-P(1)*PC4) ) 


DAV20098 


74 


A(3,3)=P42+P32-P22-P12 


DAV20099 


75 


RETURN 


DAV20100 


76 


END 



DAV20101 

r 1 

w . 

t 

i 

77 


SUBROUTINE EHPOT(P,QT) 


DAV20102 


78, 


REAL*4 P<4 ) , QT (4,4) 

43 

DAV 20103 

1 


C 

****THIS SUBROUTINE SETS UP THE MATRIX Q 

TRANSPOSE FOR ANY GIVEN 

SET DAV20104 

; 


c 

****0F EULER'S HOMOGENEOUS PARAMETERS 

the p sue I 

DAV20I05 

i 

79 


QT ( 1 t 1 ) =P ( 4) 


DAV20106 


8 0 


QT ( 1,2 )=-P (3) 


DAV20107 


81 


0T(1,3)=P(2) \ 


DAV20108 


82 


QT ( 2 r 1 ) = P ( 3 ) l 


DAV20109 

05 1 

83 


QT<2»2)=P(4) 


DAV201I0 

— — © : 

H I 

84 


QT ( 2 » 3 ) =-P ( 1 ) 


DAV20111 

W rt ! 

1 

85 


QT ( 3 » 1 ) =-P ( 2 ) * 


DAV201 12 

^ 3 

86 


QT(3»2)=P{1) 


DAV201 1 3 

O j 

• * 

97 


QT (3t3)=P(4) 


DAV20114 

KQ | 

88 


QT (4,1 )=-P( 1) 

■ 

DAV201 1 5 

O ? 

8 9 


QT ( 4 * 2 ) =-P { 2 ) 


DAV20116 

L*J ; 

** * ! 

90 


QT { 4» 3 ) =-P ( 3) 


DAV201 1 7 

Ln i 

91 


RETURN 


DAV20I 1 8 


92 


END 


DAV20119 




109 ~ SUBROUTINE NORM (P, MAG) ' " " ‘ DAV20140 

110 __ _ REAL** P { 4 ) » MAG ; DAV20141_ 

C ~ **** THIS ROUTINE NORMALIZES P TO LENGTH = UNITY AND RETURNS ITS DAV20142 

C ***»L E NGTH PRIOR TO NORMALIZATION AS A CONVENIENT ERROR MEASURE,, DAV20 143 

111 MAG = 0« DA V2 01 44 

112 DO 1 1 = 1,4 ; DAV 20145 

11 3 1 MAG=MAG+P ( f j * P ( I ) D A V2 01 46 

114 MAG= SORT (MAG) DAV20147 

115 DO 2 1=1,4 • DAV20148 

116 2 P( I )=P( I ) /MAG DAV20149 

117 RETURN DAV201 50 

118 END DAV20 151 


93 SUBROUTINE WDOT (W,WD) DAV 20 120 

C EQUATIONS OF MOTION OF OUR BODY DAV20121 

C W D ' S A R E OM'E G A DERIVATIVES? T'S ARE TORQUES, C*S AND COMM • S ARE CONSTDA V201 22 

C AND MOMENTUM BIASES D A V 2 0 123 

94 REAL*4 W(3), WD(3) DAV20124 

95 COMMON/WDOTS/ T(3), COMM(IO), Cl, C2, C3, C4, C5, C6 DAV20125 

. 96 WOTl j = T ( 1 ) *C4+C 1*W ( 2 ) *W < 3 ) -COMM ( 1 f*W ( 2 ) DAV201 26 

97 WD ( 2 ) =T (2)*C5+C2*W( 1 ) ( 3) +COMM ( 2 ) *W ( 1 ) DAV20127 

98 WD ( 3 )=T ( 3 ) *C6 + C 3 * W ( 1) * W f 2 ) “ DAV20128 

99 RETURN D AV20129 

100 END DAV20130 


SUBROUTINE MAT 1 { A , V 1 , V2 ) 
UTIL ITYMUL T IPL ICATION— USE D 
REA L*4 A (4 , 4 ) , VI ( 3) , V2(4) 
DO 1 1=1,4 
V2(I)=0o 
DO 1 J=1 , 3 

1 V 2 ( I ) = V 2 ( I ) +A ( I , J ) * VI ( J ) 
RETURN 

END “ 


RELATING OMEGA 1 


DOTS 


DAV20131 
DAV201 32_ 
DA V2 01 33 
DAV20 134 
DAV201 35 
_DAV2 0136_ 
DAV 20 137 
_DAV20138_ 
DA V2 01 39 
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FORM NO PRlNIfO 


119 

120 


SUBROUTINE RK4 " ” ' ' ' ' " ' 

COMMON /RK/ ...P { 4) , W ( 3 ) , QT ( 4, 4 ) , DT , WG 


DAV20152 

DAV20153 

— 

121 

122 


K E A L *• 4 POO (4) *0(4) — — 

REAL *4 01(3) ,02(3) , D3( 3) ,04(3) * S ( 3 > , PDH4) , PD2(4) , PD3 ( 4 ) , PD4( 4 ) 

DAV20154 

DAV20155 



C 

**** I H i S SUBRUUi INE INTEGRATES W (BODY REFERENCED ANGULAR 

VELOCITY) 

DAV20156 



C 

****AND P (THE EULFR HOMOGENEOUS PARAMETERS) FROM TIME T 

TO T+DTc 

DAV20157 



c 

****TECHNIQUE IS A FOURTH ORDER RUNGE KUTTA SCHEME. 


DAV201 58 



c 




c 

’WKtguiKfci, SUBROUTINES WDOT (USER SUPPLIED), MAT 1 , AND 

ARRAY QTo 

DAV20160 


123 


CALL EHPQT { P, OT ) 


DAV20161 


124 


CALL MAT 1 ( QT , W, PD1 ) 


DAV20162 


125 


CALL WDOT ( W,D1 ) 


DAV201 63 


126 . 


DO 10 1=1,3 


DAV20164 


127 


10 S( I )=W( I )+0o5*DT#Dl(I ) 


DAV201 65 


128 


DO 11 1=1,4 


DAV20166 


129 


•11 D(I)=P(I )+0o5*DT*PDl(I) 


DAV20167 


130 


CALL EHPQT { D , OT ) 


DAV2016R 


131 


CALL MAT1(QT,S,PD2) 


DAV201 69 


1 32 


CALL WDOT. ( S , D2 ) 


DA V 20 1 70 


133 


DO 20 1=1,3 


DAV20171 


134 


20 S ( I ) =W ( I ) +0o 5*DT*D2 ( I ) 


DAV20172 


135 


DO 21 1=1,4 


DAV20173 

i 

1 36 


21 D( I )=P ( I )+0 o 5*DT*PD2( I ) 


DAV20174 

i 

• 137 


CALL EHPQT ( D, QT ) 


DAV20175 

j. 

138 • 


CALL MAT! ( OT y $ t P03 ) 

c 

DAV20176 


139 


CALL WDOT ( S , D3 ) 


DAV20177 


140 


DO 30 1=1,3 


DAV20178 

j - 

141 


30 S( I )=W( I )+ DT*D3 ( I ) 


DAV201 79 


142 


DO 31 1=1,4 


DAV201 80 


1 43 


31 D( I )=P( I )+ DT*PD3( I ) 


DAV20181 

►d » 
pj CD 

1 44 


CALL EHPQT (D,QT) 


DAV20182 

— ' 

145 


CALL MAT1 ( QT , S , PD4 ) 


DAV20183 

CD 0 
H 

146 


CALL WDOT { S »D4 ) 


DAV20184 

td rt » 

l 

147 


DO 40 1=1,3 


DAV20185 

oi 25 

143 


40 W(I > = W( I ) + DT*(Dl( I ) + 2o*(D2<I)+D3( I) ) + D4<I))/6 0 


DAV20186 

• 

149 


PDO( 1 ) =+WO*P( 2) 


DAV20187 


1 50 


P D 0 ( 2 ) =-WO*P(l ) 


DA V 2 01 88 

Q 

_ 151 


PDO( 3)=-W0*P(4) 


DAV20189 

CO • 

152 


PD0(4)=+WU*P(3) 


DAV20190 

Ln 

***ERR0R* 

4* 

MO-O 




1 53 


DO 50 1=1,4 


DAV20191 

i 

154 


50 P< I )=P ( I ) + DT* ( PD1 { I ) +PD4( I )+2 0 *PD? ( I H-2o *PD3( I )+6o *PDO ( I ) ) /12» 0 

DAV20192 


• 155 


RETURN 


D A V 20193 

i 

156 


END • 


n 4 w->m o' 



190 

c 

subroutine' Torque (t ,b» const) 

IOEAL MAGNET SWITCHER 



DAV20233 
DAV20234 



~ C 
r 

~ SETS MAGNET POLARITIES OPPOSITE 
ur,w riF R nnTS IS DETERMINED BY 

IN SIGN TO 
COMPARING 

B DOTS 

PRESENT ( B * S ) VSo PREVIOUS 

DAV20235 
DA V 202 36 

i 


C 

c 

(BOLD'S) BODY COMPONENTS OF MAGNETIC FIELD 
T HP M l TnRQUE= M CROSS 8 


DAV20237 
DAV20238 


191 

192 


REAL*4 T ( 3 ) » B ( 3 ) » Ml, M2, M3 
B 1= B ( 1 ) 

,, BIOLD/Oo 

0/, B20LD/0o 0/ , B30LD/0o0/ 

DAV20239 
DAV20240 

c 

l 

193 
1 94 


" B 2 = B ( 2 ) 
B3 = B ( 3 ) 



DAV20241 
DAV20242 

! 

195 

, 1 QA 


BD1NEG=B10LD-B1 

BD2NEG=B20LD-82 



DAV20243 

DA V 20 246 

f 

i 

197 

1 9ft 


B03NEG=B30LD-B3 
Ml =SIGN(C0NST, BD1NEG ) 



DAV20245 

DAV20246 

j 

? 

— ! 

199 

200 


M2= SIGN (CONST, BD2NEG) 
M3=S I GN( CONST , BD3NEG) 


o 

DA V20247 
D A V 2 02 48 

i 

i 

201 
? no 


T ( l ) = M2*83-M3*B2 
T ( 2)=M3*B1-M1*B3 

> 


DA V20 249 
DAV202 50 

i 

203 

204 


T( 3 )=M1*B2-M2*B1 
B 1 OLD= B1 



DAV202 51 
DAV20252 


20 5 
20 6 


B20LD=B2 

B30LD=B3 



DAV20253 

DAV20254 

i 

i 

207 

208 


RETURN 

END 



DAV 202 55 
DAV20256 















V w 

(U (D ■ : 

157 

158 


SUBROUTINE MULVEC ( A , VI , V2 ) 
REAL*4 A (3,3), Vl(3), V 2(3) 



DAV201 95 
DAV20196 

(D 0 
H i 

CCLrt 

159 

c 

MULTIPLIES A*V1=V2 
DO 1 1=1,3 



DA V2 01 97 
DAV20198 

1 ! 

-j 2 

O > 

H •; 

160 

161 


V2( I ) =0o 0 
DO 1 K = 1 ,3 



DAV20199 

DAV20200 

L.D 

G> \ 

1 62 
163 


1 V2 ( I ) =V2 ( I )+A( I,K)*V1(K> 
RETURN 



DAV20201 

DAV20202 

OJ 

^ ; 

Ln 1 

164 


END 



DAV20203 

f 



PRINTED 


165' 

r 

SUBROUTINE BFLDSH INCLIN, A, 8) 

dptiirnr MAOMFTTO FIELD FOR VARIOUS POSITIONS IN 

ORBIT— 

SEE SEPARATE 

DAV20204 

DAV20205 


166 

c 

DOCUMENTATION 

REAL 8 { 3 ) » INCLINt PI/3ol416/ 




DAV20206 

DAV20207 

; 

167 

168 


RAD0EG=PI/180o 

BN0M=loCllE + 12/(4c.0*PI*A*A*A) 

- 



DAV20Z08 

DAV20209 


169 
1 70 


SN=0o 2026 
CN = 0, 980 




DAV202 10 
DAV2021 1 

I 

171 

172 


sI = s IN( INCL j n « rad D EG ) 
CI = CCJS ( I NCL I N^RADDEG ) 




DAV20212 

DAV20213 


173 

174 


RETURN 

ENTRY BFLDS2( EARTHRtTHETAt B) 




D AV202 14 
DAV2021 5 



c 

r, 

THETA” IS THE ANGLE IN ORBIT FORWARD 
FARTHR= W(EAP.TH)*T- LONGITUDE OF THE 

FROM ASCENDING NODfc 
ASCENDING NODE — 

(DEGREES) UAV2U2LO 

(ALSO DEGREES )DAV202 17 

1 

17 5 
176 


U=<69-, 0-EARTHR)*RAD0EG 
CU=COS (U) 




UAVtUt Jlo 

DAV20219 


177 
1 78 


SU=SIN(U) 

B(3 )=BN0M* (CI*CN+SI #SU*SN) 




DAV20220 

DAV20221 


179 

180 


C1=-CU*SN 

C2=C I *SU*SN-S I *CN 




DAV20222 
DAV20223 


181 
1 82 


10 V= THET A^RADDEG 
CV=C0S ( V) 




DAV20224 

DAV20225 

j 

183 

184 


SV = S I N ( V ) 

B ( 1 ) = 2 0*BNOM* (C1*CV+C2*SV) 

* 



DAV20226 

DAV20227 


185 

186 


B( 2)=BN0M* (C1*SV-C2*CV) 
RETURN 




DAV20228 

DAV20229 


‘ 187” 

! 188 


ENTRY BFLDS3(THETA»B) 
GO TO 10 




DAV20230 

DAV20231 


189 


END 




DAV20232 

iQ 'O 
ft) 0 

hf 








td rt 







: 

O' " ; 


txt 

o 

UJ 
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DUAL SPIN EQUILIBRIUM CONSIDERATIONS 


e 
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In this appendix we consider the equilibrium and stability of 
a dual spin configuration. Although this work has been done before 
(e.g. reference 6), we derive the results here in a simpler way. 
These results were obtained independently by the author before his 
acquaintance with reference 6 and are included here for complete- 
ness . 


Consider a body with a single momentum bias wheel; the wheel 
spins at constant rate and carries angular momentum h = {0,0,h}. 
The body is spinning about the same axis as is the rotor; for 
definiteness this is taken as the " 3 " axis. Suppose we consider 
the torque free motion of the body in the "neighborhood" of an 


equilibrium point, namely 
W 1 = w 2 = ® > w 3 = ^ 

(Cl) 

In the neighborhood of this point 

• • 

= 5 f CO 603^ 

(C 2 ,a) 

0)2 = 

6 032 / o) 2 = <So) 2 

(C 2 ,b) 

11 

m 

3 

• • 

t ($ co 3 f 0)3 - 0033 

(C 2 ,c) 

Now , 

the equations of motion/ adapted from eq. 

(1) are 

z .l“l 

+ 0)2{h + 0)3(13 - I 2 )} = ® 

(C 3 , a) 

• 

i 2 0) 2 

- o) 1 (h + 03^(13 - I^)} = 0 

(C 3 ,b) 

I 3 t0 3 

0 

11 

1 — t 
H 

1 

CM 

H 

CM 
3 
1 — i 
3 

+ 

(C 3 ,c) 


Substituting eq. (C 2 ) into these, and dropping second order terms 
(involving products of the Sio's), we obtain 

IlSwi + 6 ui 2 h + (^ + 6033) (I3 - I2).} = 0 (C 4 ,a) 

I 2 6oi 2 “ <5o)]_{h + (fl + 6003) (I3 ~ I}.)} = 0 (C 4 ,b) 

1360)3 =0 (C 4 ,c) 

' . 

Thus we see immediately 6033 = 0 , therefore 0)3 = — constant. 

Then, equations a and b can be combined to yield 

{h + fi(I 3 - I 2 )}{h + ft(I 3 - I 2 )} = 0 (C 5 ) 

Ill2 




C-l 



Report No. 90345 
Page C-2 


What we have obviously obtained here is the equation for a 
simple harmonic oscillator. Consider one special case first. If 
ft is "small", then 101 (and u >2 r which has the same form of solution) 
oscillates with circular frequency 



If ft is not small, let us define an equivalent inertia 

I 3e = z 3 + h /ft (C7) 

Then eq. (C5) reads 

StO! + 6coift 2 (I 3e “ ^ X 3e " I l ) = 0 (C8) 

Ill2 

It is clear that if either 

l3e > Ii and I 3e > I 2 (C9,a) 

or 

I 3e < X 1 and ' *3e < x 2 (C9 ,b) ^ 


that is, if I 3e is a minimum or maximum inertia, then, our system 
is stable in the Liapunov sense. On the other hand, if I 3e is an 
intermediate inertia, the system is unstable. These results are of 
a physically appealing nature since they are of the same form as 
the results for spinning rigid bodies. In fact, if in eq. (C7) 
h = 0, I 3e = I 3 , and those results are obtained directly. 

For non-rigid (i.e., damped) bodies, these requirements are 
modified somewhat. First, there is generally little damping on 
the rotor, and that which is present has a destabilizing effect 
(almost always) . Consider then a system with negligible damping 
on the rotor and a non-negligible amount on the main body. If 
the spin directions are such that the momenta of rotor and main 
body add (ft>0) , it can be shown that it is necessary and sufficient 
for stability that I 3e be the major principal axis (maximum) . Con- 
versely, if the momenta oppose each other (ft<0) , it is necessary 
and sufficient that 

h c = h + l 3 ft>0 (CIO) 

-h_ <ft<0 (Cll) 

13 

There is one additional region of stability for this system, and 
it exists only if I 3 is a major axis. The region defined by 

ft < ^h 

I 3 -max(I 1 ,l 2 ) 


(C12 ) 
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corresponds to l 3 e being a major principal axis while 9 < 0 . As 
h— ►O the above results reduce to the simple major spin axis, 
requirements . 

Definition (C7) may be combined with this observation into 
a simple rule, to wit 

h Q > IQ (CIO) 

where h Q is the system momentum bias (equal to h + 13 ft) , and I = 
max(Ii, I 2 ) • This same conclusion is obtained in reference 6 
(page A- 82 ) . 

The result (CIO) has a direct application to the attitude 
acquisition phenomena noted under the section "Terminal Orientation 
Namely, since the magnets will drive the satellite at 2wo in the 
terminal orientation, and since we wish the terminal attitude to 
correspond to a stable condition, it is necessary that (CIO) hold 
for 9 = 2too* For a given orbit (which determines wo) and given 
set of inertias, this requirement then specifies a minimum h. The 
validity of this conclusion is supported by a number of simulation 
runs. In fact, it was observed that for good terminal attitude 
acquisition performance the left side of (CIO) should be at least 
three or four times the right side. 
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APPENDIX D 


SYSTEM EQUILIBRIUM POINTS 


ft 
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Consider the equilibrium states of the acquisition system 
described in this report under two conditions. (a) The magnetic 
field is fixed in trajectory axis coordinates, that is 

df PACE (B) = 0 (D1) 

dt 

and (b) , the magnet controllers employ a linear law 

M = -KB (D2) 

• 

Where B is the time rate of change of B as measured in the vehicle. 
As mentioned in the body of this report, under the section titled 
"Stability", this system has no equilibrium points unless some 
modification is made to the flip-flop magnet control law. Now, 
let us consider the equilibrium points possible under the two con- 
ditions above. 


The behavior of our system may be completely characterized by 
the six state variables o 3 ]_ , 012 , 003 , B]_, B 2 , B 3 . Equilibrium points 
may be identified by setting the time derivatives of these variables 
to zero and solving the resulting six equations in six unknowns. 

That is , we require 

• / 

u = 0 (D3) 

B = B x to = 0 eq. (8)=(D4) 


Starting with equation (1), we make the following substitutions: 
t = Q j( torque free motion), h = 0 (constant speed wheel), and eq. ( 1 ') 

H = I* to + h. We then obtain 

1*03 = -w x (I*o) + h) = 0 (D5) 

The right hand side of (D5) equals zero according to (D3) . An 
equilibrium point must then, by definition, satisfy equations (D4) 
and (D5) . (D4) simply requires that 

o) = a(t)B (D 6 ) 


where a(t) may or may 'not be zero and may or may not be constant. 
Now, if 

h = {0, 0, h- } = h , (D7 ) 

That is, we have a pitch wheel only, then (D5) becomes, in scalar 
form - 


032033(13 - 

I 3 ) + ho32 = 0 

(D 8 , a) 

modified) 



03 1 03 3 (Ii - 

I 3 ) — ho3}_ = 0 

(D 8 ,b) 

w l w 2 ( J 2 “ 

II) = 0 

(D 8 ,c) 


D-l 
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From the above discussion we now obtain the following results 

(I) to = {0 , 0 , 0} , B = any, a = 0 is an equilibrium 
point, since it clearly satisfies eq. (D4) & (D 8 ). 

(II) o) = {0, any, h/(l 2 - I 3 )}, B || w is an equilibrium 
point. 

(Ill) o) = {any, 0 , h/dj. - I 3 )}/ B || 03 is an equilibrium 
point. 

(IV) w = {0, 0, any}, B || w is the last equilibrium point. 

These four classes contain all the equilibrium points of this system. 
In passing I note (for the mathematicians among us) that there are 
an infinity of equilibrium points. Note also that at equilibrium 
points, § = 0 so M = 0 anc ^ so t ^ ie f act that we set t = 0 at the 
outset does not affect the validity of the results obtained here. 

There is, curiously, an alternate method of obtaining these 
same results. Starting with the equation of motion (1), the control 
law (D2) , and the assumptions (Dl) and (D7) , we have 

x 1 = Il<I)]_ + io 2 h (D 9 ,a). 

T 2 = 12^2 " w l h 
x 3 = I3U3 

where we have dropped the second order terms involving products of 


the w's. Using eq. (D 2 ) and (D 4 ), we obtain 

Mi = K(o)2 B 3 - co 3 B 2 ) (D 10 , a) 

M 2 = K(o3 3 B 1 - w^B 3 ) (D 10 ,b) 

M3 = K(o3^B 2 - 0336^) 

Remembering eq. ( 5 ) ? 

x = M x B ( 5 ) 

and substituting t into eq. (D 9 ), we obtain 

-V. 

1 1001 + = B 2 K(a)3B^-— 03^63) - B 2 K(u 33 _B 2 — o^B^) (Dll, a) 

I 2 W2 “ oo^h = BiK(o 3 ]_B 2 - 03 2 Bx) - B 3 K(u 32 B 3 - U33B2) (Dll,b) 

I3W3 = B2K(w 2 B3 - 003B2) _ BiK(co 3B]^ - 10-^63) (Dll,c) 
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Using the definition of equilibrium points 

w = o 

(here_we are considering only w as a state variable, 
not B) 

we obtain three linear homogeneous equations in the three unknowns - 
w l» w 2 » and to 3 . For a solution to exist, the determinant of the 
coefficients must equal zero. That is 


-< B 2 + B 3 } 


(BiB 2 + h/K) 


B 1 B 3 


(B-^ - h/K) B]_B 3 


-< b i + 


b 2 b 3 


This equation reduces to 

“4 (B 1 + B 2 > = 0 - 


or 


B = {0, 0, Bo) 


b 2 b 3 

(B 1 + B 2 ] 


= 0 (D12) 


(D12 ' ) 
(D13) 


Substituting this back into eq. (Dll) we have 

. 2 

Iltoi + to 2 h = -B^Kto^ 

a 9 

12^2 ” = — B^Kt02 

I3OJ3 = 0 


(D14 , a) 
(D14 ,b) 
(Dl4,c) 


Since to = 0, we obtain the following equilibrium condition 

to = {0 , 0 , any) (D15) 


This result is not valid if to 3 is large, since we neglected 

<^ 3(13 “ I 2) an< 3 w 3 (I 3 - i]_) in comparison to h in writing eq. (D9). 

To the extent that this is valid, (D15) is correct. Note that (D15) 
is the same as condition (IV) identified previously, and that (I) 
is a special case of (IV). 

***>-.. 

The only special thing about this alternate derivation is that 
we did. not assume g _ q at the outset. We were considering cases 
where “ = 0 but 3 and hence t not_necessarily equal to zero. Note 
that, in retrospect, we now have B = {o, 0, B 3 }, so Bi -= b 2 = M]_ = 

M 2 = 0. The result is then that t = 0 anyway. 

In the next section we consider the stability of the equilibrium 
points enumerated here. 
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APPENDIX E 


STABILITY OF EQUILIBRIUM POINTS 
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In this appendix we investigate the stability of a particular 
equilibrium point, namely point IV of Appendix D. This stability 
investigation is carried out by linearizing the system equations 
about the equilibrium point and then forming the characteristic 
equation. A study of the roots of this equation then determines 
the stability of the system near the equilibrium point. This does 
not rule out the possibility of limit cycle like behavior but is 
nonetheless a useful technique. From the analysis that follows, 
certain conclusions may be drawn. These conclusions and their 
consequences are discussed at the end of this appendix. 


Consider perturbations about the equilibrium point 
co = { 0 , 0 , ft} 

B = {0, 0, B} 


(El) 


Then 


co 


= {Sco^, 6(02 / ft + 6003} 


(E2) 


B = { 6 B 1 , 6 B 2 , B + 6 B 3 } 

The following equations then apply 
B = B x co 
6 b-^ = ft 6 B 2 — B 60 J 2 

6b 2 = B6C0JL - fteB-L 

6B3 = 0 

where we have dropped second order terms in eq. 
a linear control law, i.e., 

M = -KB 


(8) 


(E3) 


(E3) . Now, assuming 


(E4) 


M x = K(B6co 2 - ft6 B 2 ) 
M 2 = K (ft6 Bi - B6 o5;l) 

m 3 = 0 

Then, since 

x = M x B 
We obtain 

xl = KB<ft6B;]_ - Bficoj.) 
- x 2 = KB(ft6B 2 - B6 co2 ) 


(E5) 


(5) 


(E 6 ) 


= 0 


E-l 
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Similarly, we may start with eq. (1) , substitute from (E2) 
and drop second order terms in the result. Thus we obtain the 
linearized equations of motion about the equilibrium point 

I 1 6d) 1 = t-l - 6co 2 (h + Q(l 3 - I 2 )) (E7) 

I 2 6oo 2 = t 2 + (h + ( 1 3 - !]_) ) 

I 3 Su) 3 = t 3 

These may be written in a way that makes them easier to interpret 
by using the following shorthand notation: 


l 3 e = h + I 3 (E 8 ) 

ft 

I 3e ft = h + 13 ft = ho (E9 ) 

h + ft (I 3 - 12) = ^ (I 3 e " 12) = ctl (E10) 

h + ft (I 3 - Ij_) = ft (l3 e - Ii) = a 2 (Ell) 


Using this notation and substituting for the Ti from (E 6 ) we obtain 


6<o 3 

= -KB 2 6(0 

af6 


II 1 

iY 

6u) 2 

= a 2 6(o-. - 

KB 2 5 co 2 


^2 

J 2 

SB;l 

— + 06 ( 0 ^ — 

B6oj 2 + 

6 B 2 

= B6(o^ + 

06 w 2 - 


(E12) 


The last two equations are from (E3) . Note that the last equations 
in (E3) and (E7) immediately yield 


6 B 3 = 0, 6 B 3 . = const = 0, B 3 = B (E13) 

x \ 

6003 = 0 , 6(03 = const ^^ 0 , co 3 = ft 

Equations (E12) constitute four first order ordinary linear dif- 
ferential equations in four unknowns. They are of the form 

• X = A X 

where A is a square matrix with constant coefficients. 


(E12 ,a) 
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Now, if A is negative (positive) definite it may be immediately 
shown that the system described by (El2,a) is assymptotically 
stable (unstable) by applying Liapunov's (Chetayev's) stability 
theorem to the function v = x-x. Unfortunately, for our problem 
as for many others, A is sign indefinite and the theorems noted 
are not so easily applied. 


Instead, we Laplace transform (E12,a) 



sX (s ) — x 0 = AX (s ) 

(E14 ) 

Thus 

X(s) = (A - Is ) (~X q ) 

(E15) 

and 

x(t) = L -1 { (A - Is) -1 (-Xo) } 

(E16) 

Note 

x(t) 

that_, as shown by E15, all the components of X(s) (remember, 
and X (s ) are column vectors) have a denominator D(s) and 


D(s) = Determinant (A - IS) 

(E17) 

Now, 

D(s) may also be written as 



D (s ) - (s - si) (s - S2 ) • • * (s - s n ) (E18) ■ 

where the si are the roots of D(s) = 0 or Det (A - I s) = 0 . 
Alternately, the si are the eigenvalues of the matrix A. Obviously, 
for stability, it is necessary and sufficient that the real parts 
of the s i be less than zero. 


To determine the si we expand 
D (s ) = Det (A - Is) = 0 


(E19 ) 


obtaining the (linearized system) characteristic equation in terms 
of the various system parameters. The matrix A from (E12) is 


•KB 2 /!! 

a 2/*2 


A = 


0 

B 


-ai/ 1 ! 

-kb 2 /i 2 

tB 

0 


KBft/l! 

0 


KBft/12 


fi 


-n 


(E20 ) 


an, 3 the expansion (E19) yields, after a little algebra 
Ill2S 4 + (Ii + I 2 )KB 2 s 3 + (I 1 I 2 n 2 + Ct-]_a 2 + K 2 B 4 )s2 
+ KB2Q(ai + ct2)s + aia 2 fi 2 = 0 


(E21) 



Report No. 90345 
Page e-4 


We require, for stability, that all the solutions of (E21) lie in 
the left half-plane. We may attack (E21) either numerically or 
analytically. Consider the analytical approach first. 

If we apply Routh's Criterion to (E21) we obtain, as the 
necessary and sufficient conditions for stability, that the 
following quantities must be greater than zero if aq is 


Ill2 > 0 ( al ) 

(II + 12 ) KB 2 >0 ( a2 ) 

(|)1 = Iil2ft 2 + ct]_a2 + K 2 B 4 - I 1 I 2 ^(ai + a2 ) > 0 (a3) 

II + 1 2 

4>2 = KB 2 fl(ai + a2 ) ** KB 2 ft 2 aia2 (-*-1 + ^^/^l > ® ( a 4) 

a^a2^ 2 > 0 (a5) 


Requirement (al) is trivial, and so is (a2) as long as K>0. (a5). 

is satisfied as long as the body is dynamically stable (see appen- 
dix C) without the control system. We must address ourselves to 
(a3) and (a4). We first note, en passant, that if we set K = 0 
in (E21) we obtain the characteristic equation for the system with- 
out control torques. This situation was examined in appendix C. 
Equation (E21) becomes 

s 4 + as 2 + b = 0 (E22) 

where 

b = a^a 2 ^ 2 
*1*2 

and 

a = ft 2 + aqa 2 + K 2 B 4 

Ill2 

It is trivial to show that the requirements for stability of (E22) 
are that a>0, b>0. Obviously b>0 is the stronger condition; we 
obtain the requirement 


a l a 2 > 0 (E23) 

and this is always satisfied if the system is stable without the 
control system. 
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We return now to conditions (a3) and (a4). <J>i can always be 
made greater than zero for sufficiently large gain K. For K = 0, 

(a3) may be written as 

2(11 + I 2 ) + (II + 12) U3e ~ 12) (I3e ~ H) > 2J (e24) 

Ill 2 

Various cases of this can now be considered. If A is small and 
less than zero (opposite in sense to h) , (E24) is easily satisfied. 

If ft is small and positive this relation is always satisfied. For 
large values of ft (either sign) I 3 e -®*l 3 . If I 3 is a major or minor 
axis (E24) is always satisfied. If 13 is an intermediate axis, 

(E24) may fail (for K = 0) for sufficiently large ft (which is re- 
quired to make I3e*^l3)> Before this happens the criterion (a5) 
will fail, so we need not concern ourselves with (a3) except perhaps 
for some intermediate values of ft, such that neither l 3 e >> Il^ 12 nor 
I 3 e -*l 3 . If we make some simplifying assumptions- namely Ii = 12 = 
and l 3 e = nl, (E24) reduces to 

(n - l ) 2 > (n - 2) (E25) 

which holds for all values of n, positive and negative. It would 
seem that except for a rather oddly configured satellite that <P 1 > 0 , 
and (a3) presents no problems. In fact, unless the ratio between 
Ii and I 2 is greater than 5:1, a3 must be satisfied. Additionally, 
for K>0, <f>i increases and (a3) presents even fewer problems. On the 
basis of this analysis we turn our attention to condition (a4) . 

(J >2 = KB 2 ft(ai + (X2 ) ~ KB 2 ft2 ai a 2 (I]_ + I2)/<J>1 > 0 d4) 

(j) 2 = KB 2 ft 2 [ (2I 3e -Ii-l2)-ft 2 (I 1 +I 2 ) ( I 3e _I ldl3e- I 2)/4>i] > 0 (E26) 

Now, from condition (a5) we see that the second term of (E26) is 
less than zero. Thus, if the first term is less than zero, this 
system is unstable. It is a requirement then, for stability, that 
I 3e be the major principal axis of the satellite (ie , I 3 e >max(li,l 2 ) ) • 

Now, for our system, I 3 is a minor axis, and for large values of 
ft, ISe^IS/ so the system is unstable. For ft>0, in fact, I3 e is 
always a minor axis (if 1 3 is) . The only stable values for opera- 
tion are positive values not "large". In examining this problem 
we now turn to a numerical study. 

a numerical evaluation, the following values are sub- 
in conditions. (a3)",^-(a4) and (a5) . 

35 slug- ft 2 (47.5 Kg-m 2 ) 

50 slug-ft2 (67.9 Kg-m 2 ) 

25 slug-ft 2 (33.9 Kg-m 2 ) 

0.5 Gauss 

h = 0.7 lbf-ft-sec (0.944 Nt-m-sec) 

K = 10 4 pole-cm/ (5 x 10 “ 4 Gauss/sec) = 2xlO~ 8 Nt-m/ (Tesla 2 /sec) 


For 

stituted 

J 1 = 
*2 = 
13 - 

B = 
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K is chosen by calculating the describing function for a flip-flop 
control law and a system oscillation amplitude of about wo* Note 
that K decreases as to grows and hence the system becomes less 
stable. Now, (a3) becomes 

4 ,-L = 50190 2 - 97.7a + 0.89 + K 2 B 4 > 0 (E27) 

(E27) has no real zeros even for the worst case (K = 0) , so we 
need not concern ourselves with it further. We next look at (a5) 
and determine the maximum value of a for which l 3 e is a major axis 
(recalling that this is required by ^ 2 ) 

h + I3 > max(I]_, I 2 ) (E28) 

a 

The result is a max = 0.0278 rad/sec. 

Finally, we examine <f>2 (expression a4): 

4>2 = - 53315a 4 + 5185 a 3 - (47.6 + 102 . 7 )a 2 + 1.89a (E29) 

^2 <f>l <f>l 4>1 

We are interested in determining whether <j>2>0 for .0278>a>0. This 
is required if we are to have any stable region of operation at all. 
This task is somewhat simplified by first examining <f>i. We see that 
<j>l is relatively constant for this range of a. In fact, for K = 0, 

(4>1 = 501502 - 97.7a + 0.89) (E30) 

<l>l(0) = 0.89 

( . 01) = 0.515 

4> 1 (.02) = 0.9 4 

<f> 1 (.0278) = 2.05 

and 

^lmin = 0* 4 18 at a = 0.00973 

K = 0 is surely the worst possible case, since as K increases, so 
does 4>i and thus <f>2 will increase if we are in the range 0<a<0.0278. 

Now, (E29 ) has only two real roots. One is obviously at 
a = 0, the other lies at about a = 0.05. Between these roots ^2 > 0. 
Thus we have ensured a stable region of operation for 

0.0278 > a > 0.0 (E31) 

One may convince oneself of the truth of this result by re- 
turning to equation (E21) , substituting values for the coefficients 
and then solving the resulting equation numerically for the roots 
si. It is, of course, a necessary and sufficient condition for 
stability that all the si lie in the left half plane. 
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Substituting into (E21) we obtain 

3225s 4 + 57.7s 3 + (3687ft 2 - 45ft + 1.14)s 2 + (-23.8ft 2 + 0.944ft)s 

+ (462ft 4 - 45ft 3 + 0.89ft 2 ) = 0 (E32) 

Numerical solutions for the roots of (E32) confirm the results 
(E31) obtained thru the Routhian analysis. 

All of the results obtained in this appendix, up to this point, 
are only valid in a rigorous sense, for motion in the "neighborhood" 
of the equilibrium point. Otherwise the linearizations made in 
obtaining the system equations (E12) are not allowable. Nonetheless, 
considerable digital simulation experience indicates that the results 
obtained in this section are in fact applicable to the real (non- 
linear) system. Let us consider then the significance of these re- 
sults. First of all, they indicate that the terminal attitude of the 
satellite must be such that the momentum bias points along the orbit 
normal, rather than opposite it. We remember that in the terminal 
condition the satellite has an angular velocity 2w 0 e n . If the satel- 
lite were to acquire "backwards", then the angular velocity 2 wo e n 
would correspond to ft< 0 , and we have seen that this equilibrium point 
is an unstable one. We have seen that there are no stable equilibrium 
points with very high angular velocities ; the greatest being ft = . 
0.0278 rad/sec. This means that, with the exception of possible limit 
cycle difficulties, this system should be well behaved and the control 
theory developed in the main body of this report provides for the 
positive despin of the vehicle. The author also wishes to note that 
he has observed no limit cycle behavior during extensive digital simu- 
lation of this system. 

The reader is probably asking, at this point, "How about the 
equilibrium points II and III enumerated in Appendix D" . The answer 
is that the author has not yet had time to perform a rigorous stabil- 
ity analysis for these points. However, on the basis of digital simu- 
lation it appears that they are not stable points unless 161=102 = 0 , 
which is of course the case treated here. 

Finally, we note that even the stable equilibrium states deter- 
mined here present no real despin problems. The investigations here 
dealt, for simplicity, with a fixed external magnetic field. In fact, 
even in the worst case (an Equatorial' orbit with the pitch wheel along 
the orbit normal) the Earth's field appears to cone about with a half 
angle of 11.7°. If the satellite is to remain in the equilibrium 
state, it too must cone about the orbit normal at frequency wo* The 
control torques maintaining the equilibrium state (i.e. damping out 
disturbances - which is how the coning of the B field appears to the 
system) are essentially dissipative in nature, since they oppose 
sensed velocities. The final result is then that even the small 
residual spin ft < .0278 rad/sec along the momentum bias will be elim- 
inated due to this previously unconsidered effect. 
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MAGNETOMETER NOISE MEASUREMENTS 


PURPOSE OF THE STUDY 

This study is part of an attempt to detect apparent changes 
in the earth magnetic field from a satellite during initial 
attitude acquisition. A system designed to accomplish this 
(Figure 1) consists of the Schonstedt RAM-5C-NB Magnetometer 
followed by a differentiator circuit. The output of the differ- 
entiator would be fed to a threshold detector of an attitude 
control system. This report evaluates the noise type and levels 
present in the magnetometer-differentiator combination. 



<oe Ov=l ^ kk. 


GENERAL OPERATION OF THE MAGNETOMETER-DIFFERENTIATOR 

The proposed magnetic aspect system (Figure 1) is composed 
of a probe and associated electronics (Magnetic Aspect Sensor) 
purchased from Schonstedt, and a differentiator circuit. The 
sensor provides a "d.c." output voltage proportional to the 
component of the magnetic field coaxial with the probe. The 
differentiator circuit not only differentiates, but provides 
attenuation over a range of frequencies outside of the range of 
vehicle rotation frequencies. Primarily, two circuits were used 
during the study, one with a corner frequency W c of 1 radian per 
second and the other with a corner frequency of 5 radians per 
second. The approximate equivalent circuits are shown in 
Figure 2 . 
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The mathematical representation is 

Y ( j W) = -jlOOW ; W is in radians per second 
(1+jW ) 2 

W c 

Figures 3a and 3b are the Bode plots corresponding to the circuits 
of Figures 2a and 2b. Zero db is defined equivalent to unity gain. 

TEST RESULTS 

Figure 4 is the actual circuit configuration of the differ- 
entiator with 1 radian/sec. corner frequency. Separate FET's 
were used. The circuit with a corner frequency of 5 radians/sec. 
employed the 2N3922 dual FET's. 


loM, 



subjected to each test. 

The outputs of the W c =l and W c =5 circuits were measured when 
these circuits were used in combination with magnetometer 3883. 
The two tests were not run simultaneously. The magnetometer was 
centered radially in a wooden cylinder which in turn was centered 
radially in an iron pipe. The pipe was then closed. A magnet 
external to the pipe supplied the magnetic field. The magnet 
produced a 2.3 volt magnetometer output. The output of the dif- 
ferentiators were recorded using a Brush recorder. 


August 27, 1971 
Page 3 


The peak to peak values obtained are 

8mV for the W c = 1 circuit, 

80mV for the W c = 5 circuit. 

Figures 5a and 5b are samples of the recordings from which the 
values were obtained. 

To determine the noise generated by the differentiator circuit, 
the output was measured with the input shorted. Both circuits were 
placed in a metal box and the outputs simultaneously measured using 
the Brush recorder. The peak to peak output was half a millivolt 
for the W c = 1 and 2 millivolts for the W c = 5 circuit. 

The magnetometers themselves do not have identical noise 
characteristics. Figure 6 is a Brush recording showing the noise 
output from the W c = 5 circuit fed by each of the three magnetometers. 
Each of the magnetometers was tested in the pipe with the magnet 
supplying the same field in each case. A considerable difference 
between the noise levels of the magnetometers is apparent. Dis- 
turbances, some deliberately created, in the laboratory verify that 
the magnetometers and the differentiator were functioning during 
the measurements . 

Two magnetometers, each in combination with a W c = 1 differ- 
entiator were tested simultaneously to show that spikes which were 
being observed on the oscioolgraph are not internally generated. 
Oscillographs obtained during the night hours also tend to confirm 
that the largest spikes are not internally generated. The night 
graphs are almost without the large spikes which are common in the 
data taken during normal working hours . 

Oscillographs were run for 1 hour with the W c = 1 differentiator 
and then for 1 hour with the W c = 5 differentiator. These recordings 
show that noise with a period on the order of minutes is negligible 
by comparison with the amplitude of the noise envelope. 

Root-mean-square (RMS) measurements were obtained to indicate 
the general level of the noise as a function of Wc. The measurements 
were obtained graphically from the Brush oscillographs. Three 
circuits were used in turn with magnetometer 3883. The magnetometer 
was enclosed in the pipe and the field supplied by a magnet external 
to the pipe. All three differentiators were of the same mathematical 
form, all had K = 100, and all had the same input capacitor and out- 
put resistor. With those restrictions, the noise levels may 'be 
considered to be a function of W c . The equivalent circuits for two 
of the differentiators are Figures 2a and 2b, the third (Wc = 0.5) 
is shown by Figure 7. 






Vertical Scale: lmV/small div. 
Horizontal Scale: 5mm/sec. 


Noise Output of Differentiator with W c = 1 rad/sec. in 
Combination with Magnetometer SN 3883. 

FIGURE 5A 



Vertical Scale: 5mV/small div. 
Horizontal Scale: 5mm/sec. 

Noise Output of Differentiator with W c = 5 rad/sec. in 
Combination with Magnetometer SN 3883. 

FIGURE 5B 



3882 3881 3883 

Comparison of Noise Levels from Three Magnetometers. 
Output of Differentiator with W c = 5 rad/sec. 

FIGURE 6 





































August 27, 1971 
Page 4 


\o 

-AAA/ - 


H K 

D> 


V- 


lootl. V O/A. 

vvvv~~l j — — 

t)\PFe^'cUT\KTOyt V4\TVi \U<. =. O.^ 

V NCaU^v^r- 1 

The results of the graphical analysis are given in Table 1. 


W c 

radian/sec. 


Magnetometer 3883 Output 
Volts d.c. 


Differentiator Output Noise 
Millivolts RMS 


0.5 

1.6 

0.6 

1.0 

2.3 

1.3 

5.0 

2.3 

8.1 


TABLE 1 

To obtain the RMS voltage graphically, the instantaneous value 
of the noise voltage is measured at 1 second intervals. The square 
root of the mean of the square of these measurements is the RMS 
voltage. 

MAGNETOMETER OUTPUT CHARACTERISTICS 

! 

Tests show that the d.c. output impedance of the magnetometer 
is 3.6K ohms. The 2 OK ohm resistance shown in Figure 2b actually 
includes the 3.6K ohm output impedance. 

At a frequency of W c = 0.5 radian/sec., the impedance of the 
10y farad differentiator input capacitor is 200K ohms, so there is 
neglibible loading from d.c. to W = 0.5 radian/sec. 

Discussions with the vendor indicate that no electrical isolation 
problems exist in connection with the output of the magnetometers. 

CONCLUSIONS 

To get some idea of the feasibility of using the particular 
magnetometer-differentiator combination studied, some assumptions 
of physical conditions will be made. 

For purposes of this study, a satellite will be assumed in a 
500 mile' high circular orbit. The RAM-5C magnetometer has a 
sensitivity of 4 volts per oersted. Assuming a maximum sea level 
field strength of 580 milliguass in the orbit plane, and an inverse 
cube variation of field strength with orbit radius; a voltage of~c 
1.6 volts maximum would be present at 500 miles. 

The Bode plots, Figure 3, give the voltage gain of the corre- 
sponding differentiators as a function of the frequency of a sinu- 
soidal input. For simplicity, assume the satellite to be rotating 
about an axis which is fixed relative to the earth. Assume we are 
interested in rotation frequencies from 0.5 radians/sec. down to 
the orbit rate, and assume the orbit frequency is .001 radian/sec. 
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According to the Bode plots, a magnetometer signal at the orbit 
frequency would be attenuated by a minimum of 20db. For the 500 
mile orbit, the peak differentiator output would be a maximum of 
±160mV. If the angle between the magnetometer probe and the orbit 
plane were such as to reduce the field strength by a factor of 0.1, 
then the signal would be attenuated another 20db and the peak output 
of the differentiator would be ±16mV. 

Suppose threshold detectors to be at the output of the dif- 
ferentiator. Based on the above discussion, the threshold levels 
might be set between ±16mV and ±160mV. Recall that the peak to 
peak noise values recorded on the oscillograph were 

8mV for the W c = 1 circuit, 

80mV for the W c = 5 circuit. 

The oscillographs show that the noise level measured at the 
output of the differentiator is dependent on the d.c. output voltage 
of the magnetometer. The reason is not known. 

The study does not conclusively define the limitations of this 
particular magnetometer-differentiator combination. However, based 
on the measurements , the prospects for this combination appear 
promising. 
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MAGNETOMETER NOISE MEASUREMENTS 
Addendum to Report No. 90293, Aug. 27, 1971 


Report No. 90293 gives the results of noise measurements 
for the Schonstedt RAM-5C-NB Magnetometer in combination with 
a differentiator circuit. As mentioned in that report, the 
noise level measured at the output of the differentiator is 
dependent on the d.c. output voltage of the magnetometer. 

Oscillographs of the noise at the output of the magneto- 
meter itself shows a dependence on the d.c. output voltage of 
the magnetometer. To obtain this more recent data, magneto- 
meter 3883 was enclosed in the iron pipe and a constant magnetic 
field supplied by a magnet external to the pipe. The output of 
the magnetometer was coupled to the preamp of the Brush Recorder 
through a lOpf capacitor. Samples of these recordings are pre- 
sented in the Figure. 

A close look at the oscillographs show the noise level to 
be higher for a magnetometer output voltage of -2.4V d.c. than 
for an output voltage of +2.4V d.c. Since all of the data in 
Report No. 90293 was obtained for positive output voltages, the 
following question arises. If data were taken for negative and 
positive instead of just positive output voltages, might magne- 
tometer 3881 be found to have greater p-p noise than 3883, 
instead of the other way around. To resolve this doubt, 3881 
was tested in combination with the Wc = 5 differentiator; the 
magnetometer was enclosed in the pipe and the field supplied by 
a magnet external to the pipe. 

The following table shows the p-p differentiator output 
noise as a function of magnetometer output voltage. The data 
is only approximate due to the difficulty of defining p-p noise. 
However, this data clearly shows that the noise levels are well 
below the 80mV p-p measured with 3883. 

Differentiator Output Noise vs Magnetometer 3881 DC Output 

Magnetometer Output W c =5 Differentiator Noise Output 

Volts DC . Millivolts p-p 


-2.4 


25 

-2.0 


25 

-1.0 


14 

0.0 


9 

+1.0 


20 

+2.0 


19 

+2.4 


35 


<L 
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CONCLUSIONS 

Noise levels vary from magnetometer to magnetometer. 

Noise levels also vary with field strength. Finally, and most 
important, if rates on the order of orbit rates are to be de-' 
tected, the results of the tests with the W c = 5 differentiator 
are inconclusive. Simulation will be needed to provide the 
answer. However, if a relay were installed to convert the 
W c = 5 differentiator to a W c = 1 differentiator during acqui- 
sition, detection at orbit rates should be feasible. 


Jay Levy 
ITHACO, Inc. 

735 W. Clinton Street 
Ithaca, New York 14850 
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Ill Electronics For Control of Variable Permanent Magnets 
1.0 Introduction 

In this section the circuitry is described for controlling 
the state of a variable permanent magnet as a linear 
function of an input variable. The input variable can 
be either a control voltage, (for control applications) 
or a time interval (for magnetic balancing applications) . 

The control circuitry consists of a Hall sensor that 
senses the state of the magnet, and in a feedback loop 
compares it with the input voltage. The high gain 
feedback loop forces the sensed moment to agree with 
the input voltage. Small deadbands are used so .that 
once the correct moment is achieved the magnet drive 
is automatically turned off. The circuitry and accuracy 
of this method are the subject of this section. 



The B-H curve in Figure 3 -1 illustrates the inter- 
relationship between a permanent bar magnet and an 
external field. 

In changing the permanent induction of the bar from A 
to D the path ABCD must be traversed. An external H 
field corresponding to point C must be applied and when 
removed the induction level in the bar will relax to 
point D. The continuous application of an external 
field less than that associated with point B, once 
the magnet has been cyclically magnetized, will not alter 
the permanent induction. The minor loop AB is designated 
the recoil energy and the slope the recoil slope. 
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The average slope of the recoil permeability is constant 
for a cyclically magnetized bar up to approximately 80% 
of the saturated induction level and thereafter extremely 
non-linear. 

For a bar magnet wound with a suitable coil to permit 
varying the induction level it is possible to predict 
the residual induction which will result after removing 
the magnetizing current. In the Varment system this is 
accomplished by sensing the magnetizing coil current 
and scaling it to cancel the recoil induction as 
measured by the Hall sensor. 

2.0 Control Electronics 

The Control Electronics for setting the state of the 
variable permanent magnet consists of a bipolar power 
bridge circuit, which supplies the necessary current 
with the proper polarity to the magnet coil, operated 
by one of two threshold circuits. See the Block Diagram 
(Figure 3-2) . 

Thresholding is achieved by two separate modes of operation 
designated Analog and Integrate Mode. In either mode the 
magnet state is continuously monitored by a Hall sensor 
and the current in the magnetizing coil is sensed for 
both magnitude and direction. The Analog Mode, which is 
the most straight forward, permits setting the state of 
the magnet in reference to an input analog voltage. 

The analog input. Hall output and the current correction 
voltage are summed, properly scaled, and the output of 
this summing amplifier operates the threshold circuits 
which have built in hysteresis. 

The amount of hysteresis does not affect the resolution 
but does limit the minimum incremental moment that may 
be obtained. Thus far the obtainable resolution has been 
found to be 2% of the maximum moment if the final desired 
moment is approached from the same direction each time. 

For random approaches the resolution has been found to 
be approximately 4%. Therefore in order to assure that 
the desired commanded moment will yield a differential 
output between the Hall voltage, the analog input and the 
recoil correction that is within the deadband of the 
thresholding circuit it is necessary to insert hysteresis 
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greater than 4% of the maximum capability. Otherwise 
the system could continuously switch the current on and 
off in a futile attempt to make the Hall voltage agree 
with the commanded analog voltage. 

In the integrate mode the Hall voltage is summed with a 
D.C. voltage utilizing an operational amplifier with 
several seconds lag. It is only necessary that the D.C. 
voltage exceed the threshold to operate the magnet but 
the amplitude and time lag determine the ultimate accuracy 
with which any desired moment can be attained. 

Upon commanding the system to the integrate mode and 
energizing the magnet driver the magnet state is auto- 
matically changed a fixed amount equal to the threshold 
in fixed increments of time until the maximum moment is 
reached. Should the desired moment fall between any 
two increments it would be necessary to command a 
momentary set voltage into the system that exceeds the 
threshold. When removed the magnet would then be driven 
to the desired state. To maintain this state it would 
then be necessary to de-energize the driver prior to 
the next incremental state change. 

; 

As an example consider a magnet capable of 10,000 p-cm 
initially at a 1,000 p-cm state with a desired new state 
of 7300 p-cm. Let the threshold be set at 500 p-cm and 
the integrator time constant and input set so that the 
state is changed 500 p-cm every ten seconds. The system 
would require 120 seconds to attain 7,000 p-cm and the 
momentary set command would be inserted at 126 seconds 
from the initial integrate command. The 7300 p-cm state 
would then be maintained for the next ten seconds during 
which time the driver would be de-energized. 



3.0 Variable Permanent Magnet Breadboard Controller 

A breadboard control system consisting of a power supply 
card and magnet control unit was built and trimmed to 
control a ten inch elliptical magnet made from P-6 
material. Schematics & layouts of these two cards are 
included in this report. 

The magnet was found to have a maximum moment capability 
of +/_ 6300 pole-cm. It was determined experimentally 
that utilizing the recoil current correction was 
satisfactory with only minor current limiting over 80% 
of the maximum capability of the magnet.lt was also 
determined that attempting to reduce the threshold to 
better than 4% of the available moment required careful 
gain matching and trimming. 

The X-Y plot show in Figure 3-3 show the control 
circuit and magnet being operated in the integrate 
mode starting from both the positive & negative moment 
to determine the linearity of the incremental moment 
changes. The threshold for this curve was set at -8%or 
approximately 400 p-cm increments. 



3-1 




K'AAs&ki&T cokttkoi_ block C3t^e)5=?/\KA 

=3-Z. ’ 









































































































































-VlOv 





.* b ® 









IV HORIZON SCANNER BREADBOARD 


1.0 Introduction 

2.0 Description of the X Scanner 

3.0 Operational Description 

4.0 Figures 

Appendix #1 Patent Disclosure - Horizon Sensor - Any Alti 
tude Satellite 

Appendix #2 Patent Disclosure - Horizon Sensor - Any Alti 
tude Satellite 


Appendix #3 
Appendix #4 


Report #90274 
Report #90276 



IV HORIZON SCANNER BREADBOARD 


1.0 Introduction 

As a part of the work under the subject contract, a 
demonstration breadboard of a new proprietary horizon 
scanning system developed by ITHACO was fabricated. The 
purpose of this breadboard is to illustrate the adapta- 
tion of this "X Scanner" to an existing Bendix Reaction 
Wheel. 

There are many solutions to solving the circular body 
(Earth) edge scanning problems. Complexity should be 
avoided. 

We think this has been done in the following description 
of the "X Scanner". The Patent Disclosure, Appendix 
#1 and #2, gives a detailed description of the scanner 
design and detector arrangement. 

The basis for this design is to keep the system to its 
simplest, lightest and most reliable configuration. In 
general, the fewer the number of parts, the less there 
is to get out of order. 

2.0 Description of the X Scanner 

This system consists of one reaction wheel which mounts 
the optical scanner parts. This reaction wheel mounts 
the scanning double dove prism so that there is only one 
set of bearings and one scanning optical element. The 
single window, one sealing joint, is mounted rigidly to 
the case. The lens is also rigidly mounted to the case, 
as is the bolometer (or detector system) , so that the 
optics cannot move or jump in relation to each other. 

The scan prism drive reaction wheel is oriented so that 
its axis of rotation is parallel to the flight axis, the 
center of the earth. 

The detector system, as explained in Appendix #2, consists 
of three detectors in the shape of a right triangle with 
a flake at each apex. This is required because the double 
dove prism, being set at 45° or close to it (this can vary 
due to requirements) , rotates the image at twice this 
angle. It was necessary to have two flakes contact the 
earth at the same time for Sun control or compensation. 
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3.0 Operational Description 

The scan pattern is such that one leg of the "X" crosses 
the earth with 180° rotation of the reaction wheel and the 
other leg with the next 180° rotation of the wheel. With 
the prism tilted 45° to the axis, the first 180° rotation 
crosses the earth and then the second 180° rotation has the 
prism tilted 90° from the first position. In this manner 
the common flake and the number 1 apex flake cross tangent 
to the earth on the first 180° rotation, and the common 
flake and the number 2 apex flake cross tangent to the 
earth on the next 180° rotation of the prism. 

In this manner one reaction wheel, one window, one prism, 
one lens, one detector and preamp will give a complete 
scan report on one revolution of the reaction wheel. The 
line of sight contact is tangential to the earth for the 
pair of flakes. Also, the line of sight of this optical 
system is in line with the optics, which means good 
optical performance (this is an improvement over some 
previous systems) . 

The mathematical analysis is given by Appendix #3 (Report 
#90274) and Appendix #4 (Report #90276) . 

4.0 Figure IV- 1 shows the demonstration model and the double 
dove prism. The large diameter at the left simulates an 
8" reaction wheel. This shows the prism without any 
window. . 

Figure IV- 2 shows the synchronous type scanner with a small 
plain window. The prism can be quite close to the reaction 
wheel . 

Figure IV- 3 shows the elliptical orbit scanner with a dome 
to cover the wide angle of this scan. The large distance 
between the prism and reaction wheel is required for the 
line of sight to clear the reaction wheel case for the 
lower altitudes. 


Preceding page blank 
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Figure IV-1 


Figure IV-2 


Scanner Without Window 
(to show prism) 


Synchronous Orbit Scanner 
(with flat window) 



Figure IV- 3 

Elliptical Orbit Scanner 
(with dome) 
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PATENT DISCLOSURE 

HORIZON SENSOR - ANY ALTITUDE SATELLITE 


Purpose - This horizon sensor is to operate at any altitude to 

the earth from approximately 500 N. Mi. out to greater 
than 24,000 N. Mi. within an accuracy of 0.1 degree. 



FIGURE I 


Invention - This invention describes a scanner, single head 

— consisting of a detector, a lens system and a scanning 

prism mounted on an axis of rotation shown basically 
in Figure I. This scans two paths at an angle 0 to 
each other and at 8/2 to the line of flight. The angle 
0 in the sketch is shown as 90° to illustrate. 


This accomplishes the purpose of controlling both 
Pitch and Roll with one scanner (previous systems 
used two scanners) . 


Invented by ./ 

Kennard Watson Harper April 1, 1971 
Witnessed and Understood by ' 

1) A^wv\ ^"1 \ 

2 ) ^ 


^ . 5 S? 7/' 





fv. 
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The system shown basically in Figure II shows the 
various parts: 



FIGURE II 


The optical scanner consists of a bolometer (detector) , 
a lens, a filter and substrate (cut on), a rotating 
double dove prism and a dome (or window) . 


Figure II view "A" shows the scanner looking at the 
center of the earth (the intersection of the two 
scans) . 


Invented by c . 

Kennard Watson Harper April 1, 1971/^?^,^^ 

Witnessed and Understood by ' 

V ' 


1 ) / T“ ^ '7/ 

2) *h- 
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Figure II view "B" shows the scanner in a low orbit 
(at about 500 N. Mi. above the earth) looking at the 
horizon of the earth, about 45° from the vertical. 

Figure II the view at 45° from the plan view T shows 
the profile of the double dove prism and the 45° 
look shows it swung 22.5° to look at the edge of 
the earth. 

Figure III shows the trace of rays through a typical 
double dove prism. 



The index of refraction controls the prism length 
ratio to its height. The higher the index the shorter 
the prism. The prism design does have to be propor- 
tioned for its particular application. The face angles 
can be varied to accomplish this purpose. 


Claims - 


1) One scanner head has two scan paths to control 
both pitch and roll for any altitude of orbit. 


2) One scanner drive with one axis of rotation can 
be oriented to its best position. 

3) There is one complete scan (two axis) with one 
revolution of the scanner drive, which means 
longer bearing life. 


4) The use of a double dove prism placed at a fixed 
angle of 9/2 Figure I controls angles of the scan 
pattern . 



P 
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Claims - 5) The scan pattern angle 0 Figure I can be varied by 
simply changing the swing of the reflecting hypo- 
tense from the basic 45° to the desired angle. 



Invented by 


Kennard Wa-fson ^Harper Xpril 1 , 1971 
Witnessed and Understood by 


<^ 2-4 , / 97 / 





PATENT DISCLOSURE 


HORIZON SCANNER - ANY ALTITUDE SATELLITE 


This is in addition to the disclosure of April 1, 1971. 

This additional disclosure covers the special detector 
flake arrangement which is dictated by the tilting of the 
image due to the Double Dove Prism being tilted 45° to 50°. 

When the dual flake (Sun rejection) bolometer is required 
the present two flake design will not work because the Double 
Dove Prism tilts the image 90°. The dual flakes will not cross 
the horizon simultaneously for both scans. That is, that two 
flakes are tangential to the earth. 

This means that the arrangement of the flakes will be as 
shown by Figure 4 so that one pair will be tangential to the 
horizon on one leg of the crossing and the other pair will be 
tangential to the horizon on the other leg of the crossing. 

The middle flake ("C") is the common flake. 

As the lens system flops the detector image over and 
around, Figure 5 demonstrates how this works. 

The lens both inverts and reverts the individual flakes 
and whole pattern. This is shown by the intermediate aerial 
image . 


In addition - the Double Dove Prism inverts the flakes 
and image due to the reflecting surface in the middle of the 
prism. This is shown by the last aerial image. 

As the prism in its working position is tilted 45° to 50° 
from its normal position the flakes are tilted 90° to 100°. 

Figure 6 shows the Double Dove Prism tilted 45° which tilts 
the flakes 90°. If the prism were tilted 50° the flakes would 
be tilted 100°. 

By placing these flakes as shown the first two flakes such 
as the Number 2 and common flake intersect the earth, and the 
Number 1 flake (trailing) can be blanked out as it isn't 
important. On the other half of the scan the Number 1 and 
common flakes will intersect the earth and the Number 2 (trailing) 
flake can be blanked out. The intersecting pair of flakes in 
each case is very nearly tangential to the earth. The tangential 
intersection is better at the higher altitude than at the lower 
altitudes. 




Claim - That the three flake bolometer arranged as shown 
is. required for sun rejection in the Double Dove Prism 
disclosure . 

The three flakes can be rotated into other axes with the 
same result except it is felt that the logical operation 
dictates that the blanked or unused flake be the trailing one. 


Invented by 
Kennard Watson Harper 


June 2 , 


1971 


Witnessed and Understood by 
1} 


D— -"T/ 
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Y# A. 7 y 

MEMO 

TO: File 

FROM: D. Sonnabend 

SUBJECT: X Scanner Beam Analysis 

Consider a horizon scanner of the type described by Harper.* 

The purpose of this memo is to determine the direction of the 

beam relative to the spacecraft at any position of the shaft 

B 

or wheel. Assume a coordinate system with base vectors Eg 
fixed to the case. Suppose that the direction of the telescope 
as seen from the prism is E®. Also suppose that the shaft axis 
on which the prism rotates is Eg. 

Now, for the purpose of analysis, the double dove prism acts 
like a mirror whose reflecting surface is the total internal 
reflecting surface of the prisms. Both sides of this "mirror" 
reflect. We suppose that the mirror is fixed in a rotating 
coordinate system with base vectors E™ which must obey the 
relation : 

ES=Rag(n)Eg (1) 

Within E w , the mirror is completely described by its unit 
normal U, which for definiteness we take as: 

U = E^fdX + E 3 SA = Rga (X) E^ (2) 

\ 

This definition is completely general, as any other arrange- 
ment at the same cant X differs only by some fixed phase angle 
added to n . The mirror cant X is suggested as tt/ 4 by Harper, 
but it is not necessary to pin it down in this analysis. 

*Patent Disclosure 4-6-71, by Kennard W. Harper 
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We define the beam direction by the unit vector V, and note 
that it is determined by reflecting the telescope Ej off the 
"mirror". Thus V must obey the vector relation: 

E? x U = U x V (3) 

Observe that replacing U by -U leads to the same result, 
so that both scans are included in the solution. It is easiest 
to express each vector in E . Accordingly: 

. u ?= R j!a<n> u S - *2 V U,R 2atn) 

(4) 

U B = f-cAsr), cAcri , sXj 

The general cross product relation for C = AxB is 
C a = £a$Y A f3 B Y ; and this case we have: 

e a3y^3 v y ~ e a3y^3 U y - e al3 U 3 

so that : 

^aBy v ? + E aBl )U B " 0 (5) 

Now combining (4) and (5) , a relation between the V® is 

produced for each value of a. Working this out: 

V B sX = V B cAcn (6) 

(V® + 1 ) s A + VjcAsn = 0 (7) 

(vf + l)cn + vfsn =0 (8) 

g 

The systems (.6) -(8) is homogenous in the variables + 1 , 

V 2 , and V^; so it is solvable only if the determinant of 

D 

coefficients vanishes, and then only for the ratios of the V a . 
That the determinant vanishes is easily verified; so a further 
relation between the V B will produce a solution. 
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B 

The relation E^*U = U*V is true, but is is coplanar to 
( 6 )- ( 8 ), so it is no help. 

The way out is that no use has been made of V*V = 1. 
Substituting ( 6 ) and (7) into this gives: 

(vfcAsn + sA) 2 + (vfcAcn) 2 + (vfsA) 2 = s 2 A (9) 

The solution V 3 = V® = 0 , = -1 is readily seen to be 

irrelevant. The remaining solution works out to: 

V® = 2c 2 As 2 n -1 (10) 

v B = -2c 2 Asncn (ll) 

V B = -2sAcAsn (12) 

Equations (10) -(12) are a complete description of the beam 

for an X type horizon scanner. 

It is of interest to determine the directions of motion of 

B 

the beams. To find these we find the change in V due to a 
change dn : 

dV B = 2 cAdn£cAs (2n ) , -cAc(2n), -sAcnJ (13) 

These components are the direction numbers of the instantaneous 
rotation axis of the beam. Now at either ry-=: 0:or n = it we have 
V B = £-1,0, Oj ; i.e., the beam is straight down from the telescope 
At these values: 

dAT$ (0) = 2cAdn £o , -cA , -sXj 
dV B (tt) = 2cAdrj[p^ -cA , s)J 

Since the bracketed quantities are each of unit length, the 

angle x between the scans where they cross is: 

0 2 

cx = c A -s A = c(2A) 


( 14 ) 



Thus the beam crosses itself at an angle 2X. 

As a final note, if X= tt/ 4 , the beam crosses itself at 
right angles and V B has the special form: 

V B = ~[c 2 n, sricri, sr)J (15) 

Looking down from the telescope the projection of the beam 
will look essentially like this: 





Report No, 90276 


MEMO 

TO: File 

FROM: D. Sonnabend 

SUBJECT: X Scanner Computational Analysis for Small Angles 

In a previous memo* the Harper X Scanner was analyzed 
to determine the direction of the beam relative to its outer 
case, as a function of wheel (prism) position. Here this 
analysis will be extended to find the conditions at horizon 
crossings, and use these to obtain computational equations 
for the determination of satellite attitude and altitude. 

The main results are valid only for small attitude errors. 

The notations of the previous memo will be preserved 

f 

intact, but augmented to describe satellite attitude. We 

assume a reference or trajectory coordinate system with unit 

base vectors E^ defined so that E^ is along the outward 

pointing radius vector from the earth's center, E^ is along 

the orbit normal as though the orbit were a right hand screw, 

and E^ is forward in the orbit plane to complete a right 

handed system. We suppose that the body fixed system E 

used in the previous memo is related to E by: 

E b = F „eJ (1) 

a ag p 

and for simplicity, suppose that F = I corresponds to no 
attitude error. This implies that the telescope looks straight 
down (-E^) when there is no error. 

Since a horizon scanner can give no information about yaw 
(i.e., rotation about E^) it will prove convenient to factor 
F as follows : 

F = AR 1 (ip ) (2) 

C- 

* "X Scanner Analysis " f D. Sonnabend 4-13-71 
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To assure that the rotation A contains no yaw component, we 
will require that ip be so chosen that the eigenvector of A 
is orthogonal to whenever F f I; if F = I then A = I and 
\p = 0 . 

T * 

The angle a between the beam axis V and the nadir -E^ is: 

cy = -V'eT = -V B E B ’eT 
1 1 a a 1 

and substituting from (1) this works out as: 

<»t - ~ V a F al 

Then using (2) we have: 

c Y = -v*A al (3) 

The point of using A rather than F is that the analysis is 
not cluttered with an angle which cannot be determined. 

In practice, (3) is useful whenever the beam crosses 
a horizon so that y has the special value y e given by: 

sy e = (1 + h/r e ) -1 (4) 

where h is the satellite altitude and r e is the radius of 
the earth. Various effects such as clouds and oblateness 
introduce errors through (4) in that the horizons may not 
be affected symmetrically; but these will not be considered 
in this memo. 

From the previous memo we may substitute the expressions 
for the V B in terms of the prism angle A and the wheel position 
H. When this is done each horizon crossing yields this 
relation : 

(1 - 2c 2 As 2 ti) A X 1 + 2c 2 AsncnA 21 + 2sAcAsnA 31 = cy e 


( 5 ) 
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If the attitude error is not too great, there will be 
four distinct values of n corresponding to the four separate 
horizon crossings on each turn of the wheel. Thus four, not 
necessarily independent, relations exist for computing- 
attitude. In sorting these out it is convenient to define 
the angle ri 0 by the relation: 

cy e = 1 - 2c 2 As 2 n 0 (6) 

Only the first quadrant value of n D will be used (0<n o ^/2) . 

Now for the case of no error, A = I and it follows that s^n = 
s ri 0 * The four possible solutions obviously correspond to the 
four horizon crossings, and one will occur in each quadrant. 

Thus, for no error, we can label the four n values by quadrant: 
ni = n 0 ; r \ 2 = -n - n Q ; *13 = -it + n G ; n 4 = -n 0 
The physical significance of n G is now evident — it is the 
rotation of the wheel from when the beam is vertical to when 
it crosses the horizon, if there is no attitude error. It 
depends only on satellite altitude. 

Using this null case as a guide, we can in general define 
a set of variations 6 rj ^ from null by these relations: 

hi = n o + 6 ti 1 ; n 2 = 71 ~ n o “ 5ri 2 ; n 3 = + ^o + 6r| 3 ; 

n 4 = -n Q - <$n 4 (7) 

In general, the 6 n ^ will be small if the rotation A is small. 
Returning now to A, for small rotations, any set of parameters 
will serve for computation, with equal facility. Accordingly, 
if we resolve A into pitch 9 and roll <f> by either A = R (0)R (<}>) 
or A = R 2 (<j>) R^ (9) , then for small angles we have to first order: 
A X1 = 1, A 2 i = -9, A 31 = 4> 



Report No. 90276 
Page 4 

Thus, from (5) and (6): 

4>tanAsn - esncn = s 2 n - s 2 n Q = s (n+n 0 ) s (n-n Q ) (8) 

It is clear that one of the terms on the right of (8) must 
be small. Substituting in the four definitions (7) clears 
this up, and after rejecting all second order terms the four 


equations become : 

<}>tanA - 0cn o = 26o 1 cn 0 = 2 (m - n 0 )cn 0 (9) 

<f> tanA + 0cn o = 26rucn = 2 (tt - n, - n )cn (10) 

o ^ o o 

4>tanA + 0cn o =-26n 3 cn Q = 2(n 0 - 0 3 - 7r)cn 0 XU) 

cj>tanA - 0cn o =-26n 4 cn 0 = 2 (n 4 + 0 o )cn o (12) 


The purpose of the 6ri^ in linearizing these relations 
has been served, and as they are not directly observable, they 
will now be ignored. Viewing (9-12) as a system of four 
equations in the three unknowns <p , 0 , and n Q , we see that 
there must be a relation between the observables rij_- To find 
this we compare (9) and (12) to get: 

n l “ n o = ^4 + n Q 
and also (10) and (11) : 

n ~ ^2 " n o = " n 3 " * 

This shows that r) 0 can be calculated by any of these: 

n o = - (r, l ’~Nn 4 ) = tt + l(n 3 - n 2 ) = tt + 1 (n i + n - n 4 ) 

2 2 * 2 4 1 3 4 

(13) 

Thus the internal relation between the r)j_ is: 

^1 + ^2 " n 3 “ n 4 = 2ir 


(14) 
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In practice, the amount by which (14) fails may be looked 
upon as a measure of the lack of circularity of the observed 
horizon. It is also possible to use (14) to adjust the r)j_ to 
improve the system accuracy. If 0 is the actual value on the 

X 

left of (14), then the adjusted values of the are 0 ^ (27r/rix) • 

The pitch angle 0 can be found by subtracting (9) from 
( 10 ) or ( 12 ) from ( 11 ): 

0 = tt - m - 02 = -7r “ 02 - 04 =-l (i'll + n 2 + n 3 + n 4^ 

2 

(15) 

Finally, the roll angle <J> may be obtained by adding (9) to (11) 
or ( 10 ) to ( 12 ): 

<J> = cotXcri o ( 0 i - 113 - tt) = cotXc 0 o (ir + 0 4 - 0 2 ) 

= 1 cotXc0 o (0 1 ~ 02 “ ^3 + n 4^ ( 16 ) 

In (16) the gain factor cotX is fixed by the design of the 
scanner. The factor C 0 O however, is altitude dependent. It is 
possible to calculate it from (13), or, if the orbit is 
sufficiently circular, it may be precalculated and fixed, or 
its variation can be simply ignored. To see whether the latter 
may be feasible, this table, calculated from (4) and ( 6 ) may 
be useful. The radius y e is assumed to be 6367.37km = 

3435.79nmi, and X is taken to be tt/4. 
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h-km 


200 

300 

500 




h-nmi 

0 

107.9 

161.9 

269.8 

539.6 

2698 

21,584 

Y e 

90° 

75°49 ' 

72°45 ' 

68°0 ' 

59°48 ' 

34°4 ' 

7 °54 ' 

cY e 

0 

.24491 

.29659 

.37459 

.50303 

.82839 

.99052 

^o 

90° 

60°20 ' 

57°0 ' 

52 0 16 ' 

44 °50 ' 

24°28 ' 

5 ° 35 ' 

cn 0 

0 

.49488 

.54460 

.61203 

.70925 

.91016 

.99525 


Thus, above 300km, the gain will change by less than 
two to one for any higher altitude. 






